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Abstract. We present an approximation to the Brunet-Derrida model of 
supercritical branching Brownian motion on the real line with selection of 
the N right-most particles, valid when the population size N is large. It 
consists of introducing a random space-time barrier at which particles are 
instantaneously killed in such a way that the population size stays almost 
constant over time. We prove that the suitably recentered position of this 
barrier converges at the log"^ timescale to a Levy process, which we iden- 
tify. This validates the physicists' predictions about the fluctuations in the 
Brunet-Derrida model. 

Keywords. Branching Brownian motion, Brunet-Derrida particle system. 

1 Introduction 

Consider one-dimensional supercritical branching Brownian motion (BBM): particles diffuse 
on the real line according to standard Brownian motions and split independently with rate 1 
into a random number of particles distributed according to the reproduction law (q'(A;))fcj>o, 
with mean greater than 1 and finite second moment. The physicists Brunet and Derrida have 
introduced a model of BBM with selection: Fix a (large) parameter N, and as soon as the 
number of particles exceeds N, instantaneously kill the left-most, in order to have at most 
particles at any time. This model, which we call the A^-BBM, has been studied by them and 
coauthors in extraordinary detail [20, 21, 22, 23]. In a first approximation, they model the 
system by a deterministic traveling wave of an FKPF-type equation with cutoff [20] (FKFF 
stands for Fisher, Kolmogorov, Fetrovskii and Fiskounov after [34, 46]). Assuming the validity 
of this approximation, they find that the linear speed of A^-BBM differs from the speed of 
the right-most particle in BBM without selection by a quantity of the order of log^ A^. In a 
subsequent work [22], they introduce a better, semi-deterministic approximation, which does 
not only yield a better asymptotic for the speed, but also the complete set of cumulants of 
the position of the front, all of them scaling as log~^ A^. Moreover, this approximate model 
together with numerical simulations suggests [23] that the genealogy of the system can be 
described on a timescale log'^ A^ by the celebrated Bolthausen-Sznitman coalescent [17]. 

In order to explain the presence of the Bolthausen-Sznitman coalescent, Berestycki, Beres- 
tycki and Schweinsberg [10] approximate the A^-BBM by BBM with a linear space-time barrier 
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at which particles are instantaneously killed and whose slope is exactly the conjectured speed 
of the A^-BBM (this idea was also present in [9] and indeed already in [22]). They show for 
their system that the number of particles and the genealogy of the system converge in the 
log^ N timescale to Neveu's continuous-state branching process and the Bolthausen-Sznitman 
coalescent, respectively. Here, we will push their results further and present an approximation 
of the N-BBM by BBM with a random absorbing space-time barrier (the model is defined 
properly in Section 7) with the following properties. 

1. The number of particles in the approximate model stays almost constant over the time- 
scale \og^ N . 

2. The position of the random barrier, suitably recentered, converges at the log'^ N timescale 
to a Levy process {Lt)t^o, given by 

rCO 

log Ele'^^^] = iXc + Co \ e*'^'^ - 1 - iAxl(a.<cl) A(da;). (1.1) {eq:laplace.levy} 

Jo 

Here, cq = ■y^2^^(fc — l)q{k), c e H is a constant and A(dx) is the image of the measure 
x~'^dx by the map x ^ c^^ log(l + x). 

These properties are stated in precise form in Section 7. 

In a next work, we plan to show that this approximation can be coupled with the A^-BBM 
in such a way that our results can be transferred to that model. 

1.1 Related work 

The author is aware of only two mathematically rigorous articles on the A^-BBM or the A^- 
BRW (branching random walk): Berard and Gouere [9] prove the log^ correction of the 
linear speed of A^-BRW, thereby showing the validity of the approximation by a deterministic 
traveling wave with cutoff. Durrett and Remenik [30] study the empirical distribution of A^- 
BRW and show that it converges to a system of integro-differential equations with moving 
boundary. BBM with absorption at a linear space-time barrier however is a well-studied 
process (see for example [27, 35, 37, 44, 59]) and is much more tractable than A^-BBM due 
to the greater independence between the particles and its connection with some differential 
equations [37, 54, 59]. 

In addition to its intrinsic interest, the A^-BBM is believed to be representative for general 
noisy traveling waves (see [21] or the review articles [65], Chapter 7, or [60]). There is indeed an 
exact duality relation between the FKPP equation with Gaussian white noise and a system 
of branching and coalescing Brownian motions, discovered by Shiga [64] in the context of 
stepping stone models. Recently, an estimate for the speed of a traveling wave of the noisy 
FKPP equation was established [57] which partly confirms the physicists' predictions. We 
believe that the present paper will we useful in the study of its dual branching-coalescing 
system, which could potentially lead to an improvement of the results for the noisy FKPP 
equation. 

Let us also note that branching Brownian motion without selection has a long history: 
Starting with [63] it has been studied by many authors and under various aspects, along with 
its discrete counterpart, the branching random walk. Since [55], its connection to the FKPP 
equation has raised very fruitful interactions between analysis and probability theory (see for 
example [49] and the references therein). BBM has been used in applications, for example to 
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model ecological and epidemic spread [56] or directed polymers on disordered trees [28]. During 
the last years, there has been renewed interest in the behavior of its extremal particles, be it 
the right-most only [1, 40, 2] or the whole point process formed by the particles at the right 
edge [3, 5, 6, 53]. The extremal statistics of several other models have actually been shown or 
are conjectured to belong to the same universality class as BBM, such as the Gaussian Free 
Field on a two-dimensional lattice [15, 16, 19], or the cover time of a 2D box by a random 
walk (see e.g. [29] and the references therein). 

1.2 Heuristic ideas and overview of the results 

We recall the heuristic semi-deterministic description of A^-BBM established in [22]: 

1. Most of the time, the particles are in a meta-stable state, in which the diameter of the 
cloud of particles (also called the front) is approximately L = c^^logN, the empirical 
density of the particles proportional to e~'^°^ sin(7rx/L), and the system moves at a 
linear speed "ycutoff = cq — Covr^/(2 log^ A^). This is the description provided by the cutoff 
approximation from [20]. 

2. This meta-stable state is perturbed from time to time by particles moving far to the 
right and thus spawning a big number of descendants, causing a shift of the front to the 
right after a relaxation time which is of the order of log^ N. To make this precise, we 
fix a point in the bulk, for example the barycenter of the cloud of particles, and shift 
our coordinate system such that this point becomes its origin. Playing with the initial 
conditions of the FKPP equation with cutoff, the authors of [22] found that a particle 
moving up to the point L + 6 causes a shift of the front by 

R{5) = -log (l 



CO " V L3 

for some constant C > 0. In particular, in order to have an effect on the position of the 
front, a particle has to reach Leff + 0{1), where L^s = L + Cq ^3 log log A^. 

3. Assuming that such an event, where a particle "escapes" to the point L + 5, happens 
with rate Ce"^"^ , one sees that the time it takes for a particle to come close to L^s, and 
thus causing shifts of the front, is of the order of log^ A^, which is much longer than the 
relaxation time when A^ is large. 

4. With this information, the full statistics of the position of front (the speed v and the 
cumulants of order n ^ 2) were found to be 



2 3 log log A^ 
■^cutoff ^ vr Cq- 



log-'N 

[n-th cumulant] 7r^CQ^"n!(^(n) 



t log^ N 

where denotes the Riemann zeta-function. 



(1.2) 

n ^ 2, 



Berestycki, Berestycki and Schweinsberg [10] put this description onto a rigorous foun- 
dation. They study BBM with absorption at the origin and with drift — /u, where /i = 

a/cq — vr^/Lgg = V + o{L^^). Their starting point is to introduce a second barrier at the 



-(eq: stati 
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point La = L^g — Cq a, for some large positive constant A, and divide the particles into 
two parts: One the one hand those that stay inside the interval {0,La) or get absorbed at 
0, on the other hand those that eventually hit the second barrier. This corresponds roughly 
to the division of the process into a deterministic and a stochastic part. Indeed, killing the 
particles at the right barrier prevents the number of particles to grow fast and thus permits 
to calculate expectations and variances of various quantities. The quantities one is interested 
in, for example the number of particles at the time log"^ N, will then have variances of order 
e~^, such that for large A, this system behaves almost deter ministically at the macroscopic 
scale. Moreover, the shape of the front predicted by the physicists, with a density proportional 
to e^'^"^ sin(7rx/L), follows simply from the transition density of Brownian motion with drift 
killed at the border of the interval (0, L^)- As for the particles that hit the right barrier, the 
number of descendants of such a particle will be at a later time of the order of e~^NW, where 
W is a random variable with tail P(W > x) ~ 1/x, as x goes to infinity. Moreover, the rate at 
which particles hit the right barrier is of the order of e'^/log^ A^. Putting the pieces together, 
the authors of [10] then show that the process which counts the number of particles of the 
system converges in the log^ timescale to Neveu's continuous-state branching process and 
its genealogy to the Bolthausen-Sznitman coalescent. 

In this article, we validate the physicists' predictions concerning the fluctuations of the 
position of the A^-BBM. Similar to [10], we approximate the A^-BBM by BBM with negative 
drift and absorption at a barrier, but instead of keeping the barrier fixed at the origin, we will 
make the barrier move along with process, in such a way that the number of particles stays 
almost constant. The movement of the barrier is very simple. Most of the time it does not 
move at all. Only when a particle hits a point a > and spawns a lot of descendants, we move 
the barrier to the right in order to kill particles and thus make the population size stay almost 
constant. After this system has relaxed (which takes a time of order o^), the barrier stays 
fixed again and we repeat this process, with the point a shifted by the amount the barrier has 
moved. 

Let us go into the details. Our system is defined in terms of the three positive parameters 
a, A and k, and we define A^ by a = Cg"^(log A^ + 3 log log A^ — ^4).^ Initially, we have a barrier 
located at the origin and a set of particles in the interval (0, a), such that Zq ^ kc^, where 

Zt = yae'^(^"W-'^) sin—. 

Here, we sum over all the particles u alive at time t and Xu{t) denotes the distance of the 
particle u from the barrier at time t. We then let the particles evolve as branching Brownian 
motions with branching rate 1, reproduction law q{k) and drift —fx, where fi = s/ Cq — vr^/a^. 
Furthermore, particles are killed as soon as they touch the barrier. We recall that by hypothesis 
the reproduction law q{k) has mean greater than 1 and finite second moment.^ 

The process {Ztjt^o is important for two reasons: Firstly, when a is large, the number of 
particles at a time t + s, where « s « a?, is approximately {2T:co)Ne^^Zt [10]. Hence, the 
initial condition is chosen in such a way that the number of particles is proportional to A^. 
Secondly, if we kill particles at and a, then Zt is a martingale and therefore very easy to 
handle. 

^We use the letter a instead of La for typographical reasons. 

^This last condition is only technical and we believe our results to be true for more general reproduction 
laws. In fact, in Section 4, we show that the asymptotic results on the random variable W obtained in [10] 
still hold if Yuk>i ^log^ kq{k) is finite, but we don't know whether this condition is sufficient for Theorems 7.2 
and 7.3. 
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When a particle hits the right barrier at the time T, say, we absorb its descendants at the 
space-time hne J£ = a — y + {cq — fi){t — T), where y is a large constant depending on A only 
(this idea comes from [10]). In doing so, the number of particles absorbed at the barrier has the 
same law as in BBM with absorption at a critical line, i.e. a line with slope cq. Defining then 
Z' as Zt, but summing only over these descendants, we know that at a later time the number 
of descendants of this particle will be of the order of e^^NZ'. Consequently, we say that a 
breakout occurs, whenever Z' > ee^, where £ will be chosen such that in particular £ « 1/A. 
Looking at the definition of Zt, it is easy to guess by which amount A we have to move the 
barrier in order to counterbalance the breakout: Choosing A = Cq ^ log(l + Z'/{Ke^)), the 
value of Zt is approximately divided by 1 + Z'/{Ke^), such that after the relaxation time, the 
value of Zt and the number of particles should again be approximately kc^ and {2'kcqk)N . 
This is basically true, but we also have to take into account the fluctuations of Zt between the 
times and T, which are mostly due to the particles hitting the point a without producing 
a breakout. For this reason, the actual definition of A in (7.1) differs from the one given 
here. Nevertheless, the above considerations already explain the convergence of the barrier 
to the Levy process given by (1.1): One the one hand, we have Z' « {'k/cq)W, where W is 
the random variable mentioned above, such that the law of e~^Z' conditioned on Z' > ee^ 
is approximately dx for large A and a.^ On the other hand, we will show that 

breakouts occur at a rate proportional to e~^a~^. Together with the definition of A, this 
explains the Levy measure A(dx) of (1.1). One easily checks that the cumulants of this Levy 
process coincide with (1.2). 

We want to stress two more points. First, in [10], the authors cut the interval [0, a^] into 
tiny pieces of size 0a^, with 9 « e^^, in order to make sure that with high probability at 
most one breakout occurs during a single piece. In adapting this approach to our system with 
the moving barrier, we found it however difficult to control the fluctuations of the process Zt 
over the whole interval of time [0, a^]. We therefore chose another approach, which also has 
the advantage of giving more information about the history of the particle that causes the 
breakout. Namely, we will classify the particles into tiers, according to the number of times 
they have hit the point a and come back to the space-time line ^ mentioned above. Thus, 
when a tier particle hits the point a, it advances to tier 1, and its descendants have a second 
chance to break out after having come back to We can then define the time T of the first 
breakout and will indeed show that T is approximately exponentially distributed with rate 
proportional to e^^a^^. Interestingly, we will see that with high probability breakouts only 
occur from particles of tier or 1, the number of breakouts occuring from particles of tier 1 
between the times and being approximately proportional to A (and the remaining ^ 
breakouts occurring from particles of tier 0). In order to study the system up to the time 
T, we will then study BBM conditioned to break out at time t for every t ^ 0, which can be 
formulated in terms of a Doob transform of the process. 

The second point concerns the shape of the barrier. If we were only interested in the state 
of the system after is has relaxed, we could shift the barrier instantaneously by A. However, 
since in a second work we plan to couple the model of this article with the iV-BBM, we need 
to move the barrier continuously and over the timescale a^, which is the relaxation time of the 
system. In this article, we will allow the shape of the barrier to be given by an arbitrary family 
of "barrier functions" (/a)a>Oj which are non-decreasing, twice differentiable functions with 
/a(0) = and /a(+cC)) = plus some uniformity conditions. First-moment calculations 

■^The statement "for large A and a" means that we let first a, then A go to infinity, see Section 6.1. 
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then suggest that the right barrier function to choose for the couphng with the A^-BBM is 



where 6{x,t) is defined in (2.2). 

We finally remark that although parts of this article (mostly in Sections 4 and 5) are heavily 
based on [10], it is entirely self-contained. This means that we will reprove some results of 
[10], often because we need stronger or different versions, but also sometimes because we found 
simpler proofs. We think that this is for the benefit of the reader. 
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2 Brownian motion in an interval 

{sec ibm} 

In this section, we recall some explicit formulas concerning real-valued Brownian motion killed 
upon exiting an interval. These formulae naturally involve Jacobi theta functions, since these 
are fundamental solutions of the heat equation with periodic boundary conditions. We will 
therefore first review their definition and some of their properties. 



2.1 Jacobi theta functions 

In [31], p355, the Jacobi theta function 9-^(v\t) is defined for e C, r e C with Imr > as 

Gsivlr) = ^ exp (j,7r{n'^T + 2nv)^ = 1 + 2 ^ e*''"^^ cos(27r?Z'u). (2.1) ^eqithetas} 

neZ n=l 

For our purposes, the following definition will be handier: For x e C, t e C with Het > 0, we 
define 



e{x,t) = ^3(1 ^) = 2 ( - 4^^* + ^^^^) 

QQ 2 yZ.Z J {eqitheta} 

= 1 + 2 e"^"'* cos(7rnx). 



n=l 



The definition (2.2) is a representation of as a Fourier series, which is particularly well 
suited for large t, but which does not reveal its behaviour as Ret 0. This is where the 
following representation comes in, which is related to (2.2) by the Poisson summation formula 
(see [8], §9): 

2 -i=»p(-(^L^). {2,3, 

neZ * 

One recognizes immediately that for real x and t, 9{x,t) is the probability density at time 
t of Brownian motion on the circle R/2Z started at 0. In other words, 6(x, t) is the unique 
solution to the PDE 

|u(x,t) = i(|:)'n(x,t) (PDE) 
u{x,t) = u{x + 2,t) (BC) 
.«(^,0+)=I]„,z5(x-2n) (IC), 

where S{x) denotes the Dirac Delta-function. This is the heat equation with periodic boundary 
condition and the Dirac comb as initial condition. Note that (PDE) and (BC) also follow 
directly from (2.2). 



{eq: theta.gaussian} 



2.2 Brownian motion killed upon exiting an interval 

Various quantities of Brownian motion killed upon exiting an interval can be expressed by 
theta functions. For x e R, let be the law of Brownian motion started at x, let {Xt)t^o 
be the canonical process and let Hy = mf{t ^ : Xt = y} . For a > and x e [0, a], denote by 
^killed a Brownian motion started at x and killed upon leaving the interval (0,o). 

Let Pt{x,y) be its transition density, i.e. 

p?(x, y) = W^aiedA^t e dy) = W^'{Xt e dy, Hq a Ha < t)/dy, x,ye [0, a]. (2.4) {eq:c,ef.p> 
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Then Pt{x, y) is the fundamental sokition to the heat equation (PDE) with boundary condition 

u(0,t) = u{a,t) =0, t ^ 0. 
Hence (see also [41], Problem 1.7.8 or [18], formula 1.1.15.8), 

rf(.,,)^o-'(.f^L^,-|)_,(2L±l,i^)y (2.5) 

Equation (2.2) then yields 



a a" 



a 



{eq:p_theta} 



Ptix,y) = - e *sin(7mf)sin(7rn^). 



This representation is particularly useful for large t: Define 



2 -7r2/2(n2-l)t 



n e 



n=2 



By (2.6) and the inequality | sinnx| ^ nsinx, x e [0, vr], one sees that 



2 ^^t, 



Pt{x,y) = -sm{TTx/a)sm{7Ty/a)e ^ {1 + 0{l)Et/^2). 



(2.6) 



(2.7) 



(2.^ 



{eq:p_sin> 



-Ceq:def_E} 



Note that the potential kernel is given by 

-HoAHa 



nao r-tlQAtla > 

p^(x,2/)dt = t^"(| l(x,eds,)dij/dy = 2a-HxAy)(a-xvy), (2.9) 

by the formula for the Green function of Brownian motion (see e.g. [43], Lemma 20.10, p379). 
Set H = Hq a Ha and define 



-[eq:p_estimate} 



-[eq:p_potential} 



rf{x) = W''{He dt, Xh = a)/dt. 



Then (see [18], formula 1.3.0.6), 



(2.10) 



(2.11) 



{eq: def _r} 



■[eq:r_theta} 



where 9' denotes the derivative of 9 with respect to x. 

The following two integrals are going to appear several times throughout the article, which 
is why we give some useful estimates here. For a measurable subset S R, define 



P{x,S) = W^(e^'''^l^Ho>H.,s) 



e2^^r"(x)ds, 



and 

r{x,y,S) = f 
which satisfy the scaling relations 



Sn(0,oo) 



(2.12) -Ceqila} 

(2.13) 



r (x, S) = I{x/a, S/a^), J^ix, y, S) = aJ{x/a, y/a, S/a^), 



(2.14) 



-teq: Ja> 



{eq: Unsealing} 



with 1 = 1^ and J = . The following lemma provides estimates on I{x,S) and J{x,y,S). 
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{lem: I.estimates} 

Lemma 2.1. There exists a universal constant C , such that for every x e [0, 1] and every 
measurable S cz R_|_, we have 



\I{x,S) — 7rA(S') sm(7rx)| ^ C|^l a £'infs(l a A(5)) sm(7r3;) j , and 
\J{x,y,S) - 2A(S') sin(7rx) sin(7ry)| ^ c(^[{x a y)(l - (x v y))] a E^inf s sin(7rx) sin(7ry) 

where X{S) denotes the Lebesgue measure of S and E[^fs is defined in (2.7). 

Proof. First note that I{x, •) is a positive measure on R_|- for every x e [0, 1], such that 

0^/(x,5n[0,l]) «:/(x,[0,l]) ^t^^(e^^il(H,^i)) ^ , 

by (2.12). Furthermore, decomposing I[x,S) into 

I{x, S) = I{x, S n [0, 1]) + I{x, S n (1, oo)), 

it is enough to prove that |/(x, S") — 7rA(S') sin(7rx)| ^ C(l a A(5))£^inf s sin(7rx) for all S. Now, 
by (2.11) and (2.2), 

I{x,S) = - f e^'9'{x,s)ds 
2 JS 

= TT 2 e-'^("'-^)'(-l)"-^nsin(7rnx)ds 

rt=l 

00 p 2 

= 7rA(S) sin(7rx) + ^ 2 J V("'-i)Ms)n(-l)"-i sin(^nx). 



where the exchange of integral and sum is justified by the uniform convergence of the sum for 
s ^ 1. We now have for n ^ 2, 



'^00 o O 9 



Js JinlS TT^in^-l) 



as well as 

^ . . „ _ 2 

Furthermore, we have for n ^ 2 



JS 

|n(— 1)"^^ sin(7rnx)| ^ n^sin(7rx) ^ 2(n^ — l)sin(7rx). 

It follows that 

4 

\I{x,S) — 7rA(S') sin(7rx)| ^ (— a 7rA(5'))i?inf 5 sin(7rx). 

TT 

This proves the statement about /. The proof of the statement about J is similar, drawing 
on (2.6) instead and on the following estimate: 

ri ^2 ^2 TOO ^2 

J(x,7/, [0, 1]) = e~^pt[x,y)dt ^ e~ \ pt{x,y) dt = e~ {x a y){l - {x v y)), 
Jo Jo 

by (2.9). □ 



2.3 The Brownian taboo process 

The Markov process on (0, a) with infinitesimal generator 



{sec :taboo> 



d 

dx 



vr vrx d 
— cot — 

a a dx 



is called the Brownian taboo process on (0, a). It is a diffusion with scale function s{x) and 
speed measure m(dx), where 



vr TTX 

s(x) = — cot — 
a a 



and 



m[axj = — ^ sm I — 1 ax. 



The singular points and a are therefore entrance-not-exit. For x e [0, a] we denote the law 
of the Brownian taboo process on (0, a) started from x by W^^^^^^ ^. Often we will drop the a 
if its value is clear from the context. 

The name of this process was coined by F. Knight [45] who showed that it can be interpreted 
as Brownian motion conditioned to stay inside the interval (0, a) (hence, and a are taboo 
states). When a = vr, the Brownian taboo process is also known as the three-dimensional 
Legendre process, because of its relation to Brownian motion on the 3-sphere (see [41], p270). 
Readers familiar with the 3-dimensional Bessel process will notice that it can be obtained from 
the Brownian taboo process as the limit in law when a — > oo. Note that the normalisation 
of the scale function and speed measure from the last paragraph was chosen in such a way 
that they converge, respectively, to the scale function and speed measure of the 3-dimensional 
Bessel process, as a ^ oo. 

Below we list some useful properties of the Brownian taboo process: 

1. It satisfies the following scaling relation: If Xt is a Brownian taboo process on (0,1), 
then aX^^^2 is a Brownian taboo process on (0, a). 

2. It is the Doob transform of Brownian motion killed at and a, with respect to the 
space-time harmonic function h{x,t) = sin(7rx/a) exp(7r^t/(2a^)). In other words, for 
X e (0,a), VF^aboo obtained from WjfjQg^ by a Cameron-Martin-Girsanov change of 
measure with the martingale 

/ . 7rx\-i . irXt TT^ 

Zt = I sm — ) sm exp — ^t. 

\ a / a 2a^ 

3. As a consequence, its transition probabilities are given by 

taboo(o,a). ^ A \ sin(7ry/a) N 



sin(7rx/a) 



(2.15) 



Equation (2.8) now implies that 

taboo(0,a) ^ ^ 2 



Pt 



x,y) = - sin (7ry/a)(l + 0{l)Et/^2), for all x,y e [0,a], 



(2.16) 



4. As can be seen from above or directly, it admits the stationary probability measure 

{m{0,a))^^m{dx) = 2/asin^(7rx/a) dx. 



{eq: taboo_transition 



{eq:ptaboo_estimate} 
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5. It is self-dual in the sense that for a measurable functional F and t > 0, we have 



taboo 



[F((X,;0 ^ s ^ t))] = ^fi^jF((Xt_,;0 ^ s ^ t))]. 



Here W^aboo denotes the taboo bridge from x to y of length t. This follows from the 
self-duality of killed Brownian motion. 
The following lemma will be needed in Sections 6 and 7. 

Lemma 2.2. Define k{x) = e^'^^. There exists a constant C , depending only on c, such that 
we have for every x,y e [0, a], 



{lem : k.integral} 



^^taboo 



f k{Xs) 
Jo 



ds 



^ C[t/a'^ + err(x) 



^2.1 7) {eq: taboo_integral_e 



and for t ^ , 



K'b'o'o 



J k{Xs) ds ^ C[t/a-^ + err(x) + err(y)^ , 



with err(2;) = {1 a z "'^) + (1 + z)e . If t ^ , we still have for x,y ^ a/2, 



taboo 



(\{Xs) 

JO 



ds 



^ C. 



(2.18) 



(2.19) 



Proof. We first show that (2.17) implies (2.18). By the self-duality of the taboo bridge, we 
have 

-4/2 rt/2 



^^taboo 



f k{Xs 

JO 



ds 



^^taboo 



rt/z 

k{Xs 

JO 



ds 



^^taboo 



r 

k[X,) ds 

JO 



It therefore remains to prove that 



rt/2 

Jo ' 



k{Xs)ds ^C{t/a^ + exx{x)). 



Conditioning on a{Xs]Q ^ t ^ V^)) this integral equals 

~p\f^°°[Xt,^,y) rt/2 



E{x,y) = W,l^,, 



p^rf^^°°{x,y) 



r 

k{Xs) ds 

JO 



By (2.16), there exists a universal constant C, such that for t ^ a^, 

'-t/2 



E{x,y)^CW^ 



taboo 



rt/z 

k{Xs) 

JO 



ds 



Equation (2.17) therefore implies (2.18). 

Heuristically, one can estimate the left side of (2.17) in the following way: Since k{x) is 
decreasing very fast, only the times at which Xg is of order 1 contribute to the integral. When 
started from the stationary distribution, the process takes a time of order to reach a point 
at distance 0(1) from [50] and it stays there for a time of order 1, hence the integral is of 
order t/a^. When started from the point x, an additional error is added, which is of order 1, 
when X is at distance of order 1 away from 0. Adding both terms gives the bound appearing 
in the statement of the lemma. 



{eq: taboo_integral_e 



{eq : taboo_integral_e 
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The exact calculations are most easily performed in the following way. Let y be a random 
variable with values in (0,a) distributed according to m{dx) := 2/a sin^(7r2;/a) dx, which is 
the stationary probability measure of the taboo process. Let Hy = inf{t > : Xg = Y}. We 
then have 



taboo 



f k{Xs) 

JO 



ds 



^^taboo 



f kiXs) ds + f k{X,) d 
Jo JHy 



^ W; 



taboo 
h + h 



(•Hy 

k{Xs 

Jo 



ds 



Y 

W; 



taboo 



f HXs 

Jo 



ds 



The second term is simply equal to 

ra rco 

l2 = t\ m{dy) k{y) dy ^ 2TTH/a^ 6-^^(1 + y)y'^ dy CT/a^, 
Jo Jo 

for some constant C depending only on c. 

The first term is equal to (see e.g. [61] Chapter 3, Corollary 3.8) 

nx ra na ry 

h = fh^dy) m{dz)Gy^aix,z)k{z) + m{dy) ■m{dz)Go,yix, z)k{z) 

Jo Jy Jx Jo 

= : 111 + /i2, 

where the Green functions are defined by 

Gy^a[x, z) = s(x A z) — s{y) and Go^y{x, z) = s{y) — s{x v z). 
By Fubini's theorem, the first term in (2.20) is easily bounded by 

111 sS m{dz)k{z) \ m{dy)[s{z) - s{y)], 
Jo Jo 

and noticing that sign(s(2;)) = -l(2<a/2) + l(^>a/2): we get 



hi ^ 



J^a rz ra rz 

m{dz)s{z)k{z) m{dy) + m{dz)k{z) m{dy){—s{y)) 
a/2 Jo Jo Jo 



(2.20) {eq : proof _k_integral 



^C/a3( r z^k{z)dz+ \ z^k{z)dz 

^ Ja/2 Jo 

^ C/a\ 

where again we made use of the inequality sinx ^ x for x e [0, vr]. 

For the term I12 a little bit more care is needed. Using the fact that J^m(d?/) ^ 1, we 
have 



h2 ^ 



J -a ry rx ra/2 

rh{dy)s{y) m{dz)k{z) + {—s{x) v 0) m{dz)k{z) + m{dz)\s{z)\k{z) 
X Jo Jo Jx 



= '■ lui + 1x22 + -^123- 

To estimate the first two terms, note that 



m{dz)k{z) ^ C(l A y^), and m{dy)s{y) ^ C/a. 

Jx 
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such that 



/i2i + /i22 ^ C(l/a + (1 A x^){-six) V 0)) s: c(l/a + (1 a x'^)^ , 
because —s{x) ^ 1/x for x e [0, a]. The third term is seen to be bounded by 

J -00 
zk{z) dz ^ C7(l + x)e-'=^. 
X 



Altogether, we get 



taboo 



J k{Xs) ds ^ C(^t/a^ + l/o + err(x)) , 



and the 1/a term can be dropped, because t ^ by hypothesis. This proves (2.18) and 
therefore (2.17). 

When t ^ a different method of proof is needed. First we note that for < x,y ^ a/2, 
the transition density of the taboo bridge can be written 



taboo 



Xg e dz 



p1{x,z)p'i_,{z,y) 
Pt{x,y) 



dz. 



If we denote by p^{x,y) = (27rt)^^/^ exp( — (z^ + a;^)/2t)2 sinh(zx/t) the transition density 
of Brownian motion killed at 0, then we have the trivial inequality p1{x,y) ^ Pt{x,y) and 
furthermore by scaling we see that Pf{x, y) ^ Cp^{x, y), since x,y ^ a/2 and t ^ a^. It follows 
that 



^^taboo 



f k{Xs) 

JO 



ds 



f k{Xs) 

JO 



ds 



where R^'^'f denotes the law of the Bessel bridge of dimension 3. This Bessel bridge is the 
Doob transform of the Bessel process started at x with respect to the space-time harmonic 
function hy{z, s) = p^_g{z, y)/p^{x, y). By the standard theory of Doob transforms, this is the 
Bessel process with additional drift 



d; 



-(log hy{z,s)) 



z d T zy 
h — log smh 

t — s dz t — s 



coth 



Now, this in an increasing function in y, and standard comparison theorems for diffusions (see 
e.g. [61], Theorem IX. 3. 7) now yield that for yi ^ ?/2) we have 

R^^i<y2Yk{Xs)] ^ i?^''*'^i[fc(X,)], 

since is a decreasing function. This is true in particular for yi = 0. Using the self-duality of 
the Bessel bridge, we can repeat the same reasoning with x. We thus have altogether 



taboo 



\\{X, 
Jo 



ds 



f k{X, 
Jo 



ds 



for any x,y ^ a/2. This calculation can be done explicitly and yields (2.19). 



□ 
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3 Preliminaries on branching Markov processes 

{sec :preliminaries} 

In this section we recall some known results about branching Brownian motion and branching 
Markov processes in general. 

3.1 Definition and notation 

{sec :preliminaries_d 

Branching Brownian motion can be formally defined using Neveu's marked trees [58] as in [26] 
and [25]. We will follow this path here, but with slight differences, because we will need to 
consider more general branching Markov processes and the definition of branching Brownian 
motion in [25] formally relied on the translational invariance of Brownian motion. 
We first define the space of Ulam-Harris labels, or individuals, 

where we use the notation N* = {1,2,3,...} and N = {0} u N*. This space is endowed with 
the ordering relations < and < defined by 

u < V 3w e U : V = uw and u < v u < v and u v. 

A tree is by definition a subset t f7, such that 0et, fetifu<'U and net and for every 
u there is a number e N, such that for all j e N*, we have uj e t if and only if j ^ k^- 
Thus, ku is the number of children of the individual u. We denote the space of trees by 3" 
and endow it with the sigma- field generated by the subsets J3^ = {te^:net}. 
For a tree t e ^ and u G i, we define the subtree rooted at u by 

= e [/ : TO e t}. 

Given a measurable space ^ , a marked tree (with space of marks ^) is a pair 

f^ = (t, (r?„;net)), 

where t e =^ and 7]^ £ ^ for all u e t. The space of marked trees is denoted by 
and is endowed with the sigma-field j?/-'^ = tt^^{£/), where vr : S^-'^ ^ \s the canonical 
projection. Accordingly, we also define = 7r^^(=5^). The definition of a subtree extends 
as well to marked trees: For n e t, we define 

(t^)W = (tH,(,^,,;^;etW)). 

For our purposes, the space of marks ^ is always going to be a function space, namely, for a 
Polish space £' and a cemetary symbol A ^ (f, we define the Skorokhod space D{S') of functions 
H : [0, od) <§ kj {A} which are right-continuous with left limits, with H(0) 7^ A and for which 
= A implies H(s) = A for all s ^ t. Then we define C('^) = inf{t ^ : = A}. For an 
individual u e U, its mark is denoted by H„ and we define Cu = Ci'^u)- The branching Markov 
process will then be defined on the space (we suppress the superscript D{S)) 

= = (t, {^u]ue i)) e Sr'^^^^ ■.\jueU \ll^i^ku:Cu<^^ H«(C«-) = S»(0)}, 
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endowed with the sigma-field ^ = $7 n s^^^^> generated by the sets 0^ = n ^ . We 
define for u e U the random variables 



which are the birth and death times of the individual n, respectively. We then define the set 
of individuals alive at time t by 

J^{t) = {uei:bu^t < du}. 
The position of u at time t is defined for n e t by 



Xu{t) 



Ey{t — if f e ^{t) and v < u 

A, if du ^ t. 



Now suppose we are given a defective strong Markov process X = {Xt)t^Q on S", with 
paths in D{S). The law of X started m. x e S will be denoted by P . For simplicity, we 
will assume that for every x e , we have C(-^) < co, P -almost surely. Furthermore, let 
{{q{x,k))ke'N)xe^ be a family of probability measures on N, measurable with respect to x. 
Then we define the branching Markov process with particle motion X and reproduction law 
q as the (unique) family of probability measures (P^)xe^ on which satisfies 

P''{dco)=P"{dX0)q{X0{C0-),k0)YlP^0^^0-\dJ'^). (3.1) M:b.p„i> 

Note that by looking at the space-time process {Xt,t)t^o, we can (and will) extend this 
definition to the time-inhomogeneous case. 



3.2 Stopping lines 

{sec : stopping. lines} 

The analogon to stopping times for branching Markov processes are (optional) stopping lines, 
for which several definitions exist. For branching Brownian motion, they have first been 
defined by Chauvin [25]. The definition we are giving below is equivalent to the definition 
there, although there are formal differences. Note that Jagers [42] has given a more general 
definition of stopping lines for discrete-time branching processes, and our definition of stopping 
lines is partly inspired by the exposition there. Note also that Biggins and Kyprianou [11] 
build up on Jagers' definition of stopping lines and define the subclasses of simple and very 
simple stopping lines (again for discrete-time processes). Chauvin's definition (and therefore 
ours as well) then corresponds to the class of very simple stopping lines. 

We first define a (random) line to be a set i = £{uj) c [/ x [0, od), such that 

1. ue ^{t) for all (n,t) e £, and 

2. {u,t) 6 i implies {v, s) f I for s\\ v < u and s < t. 

Note that a line is at most a countable set. For a pair {u, t) e U x [0, oo) and a line i, we write 
i < (n, t) if there exists {v, s) e i, such that v < u and s ^ t. For a subset A ^ U x [0, go), we 
write £ < A if I < {u,t) for all (u, t) e A. If ii and £2 are two lines, we define the line ii a £2 
to be the maximal line (with respect to <), which is smaller than both lines. 
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We now define for eacli ueU two filtrations on by 

^u{t) = (f]„n(T(H„(s);0^s^t-6„)) v Y(!^„na(H„)) 

v<u 

^r{t) = (a, n (7(H„(s);0 ^ s^t- 6„)) v \/ {VL, n (7(H,)). 

Informally, ^u{i) contains the information on the path from u to the root between the times 
and t, and ^S^'^{t) contains this information and the one concerning the descendants of 
u after the time t. In particular, we have J^u{t) ^S^^{t) The filtration ^u{t) is denoted 
by ^u{t) in Chauvin's paper [25], and c^S^^{t) corresponds to the pre-{u,t)-sigma-algebra as 
defined by Jagers [42]. 

We can now define a stopping line ^ to be a random line with the additional property 

3. V(n, t) e U X [0,ao) : {u e nu : ^ < {u, t)} e 

The sigma-algebra of the past of .if is defined to be the set of events E e such that 
for all {u,t) e [/ X [0,oo), 

For example, for any t ^ 0, the set ^{t) x [t] is a stopping line. If T = T{X) is a stopping 
time for the strong Markov process X, then 

= {{u,t) e [/ X [0,oo) : ue ^{t) and t = T(X„)}} 

is a stopping line as well. 

The first important property of stopping lines is the strong branching property. In order 
to state it, we define for t ^ 0, u e ^{t), 

with H^(-) = 'Eu{- + t — bu) and = for v e t^*'^\{0}. The strong branching property 
([25], [42]) then states that for every stopping line conditioned on the subtrees 
for {u,t) e .if, are independent with respective distributions p^"(*). 

3.3 Many-to-few lemmas and spines 

{sec : spine} 

Another important tool in the theory of branching processes is the so-called Many-to-one 
lemma, and its recently published extension, the Many-to-few lemma [38] along with the 
spine decomposition technique which comes along with it and has its origins in [52]. Here 
we state stopping line versions of these lemmas, which to the knowledge of the author have 
not yet been stated in this generality in the literature, although they belong to the common 
folklore. We will therefore only sketch how they can be derived from the existing literature. 

We assume for simplicity that the strong Markov process X admits a representation as 
a conservative strong Markov process X with paths in D{S), which is killed at a rate R{x), 
where i? : ^ [0, oo) is measurable. The law of X started at x is denoted by and the 
time of killing by C. Given a stopping time T for X, we can then define a stopping time T 
for X by setting T = T, ii T < and T = oo otherwise. For simplicity, we write Jfp for 
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Finally, for every x e (f , define m{x) = Xifc>o(^ ~ ^)Qi^jk), mi{x) = Y^k^^kq^x^k) and 
m2{x) = Y^k^o k{k - l)q{x, k). 

We are now going to present the spine decomposition technique, following [36]. They 
assume that g(x,0) = 0, but this restriction is actually not necessary, as noted in [38]. Given 
a tree t, a spine of t is formally an element of the boundary of t, i.e. it is a line of descent 
i = (Co = 0) Ci) '^21 • • •) from the tree, which is finite if and only if the last element is a leaf of 
the tree. We augment our space O to the space 0* by 



= {{uj, ^) : (x! e $7, ^ is a spine of the tree underlying oo] 



We are going to denote by the individual u e U that satisfies u e ^{t) and li e ^, if it 
exists^, and = otherwise. Instead of writing X^^(t), we are going to write for short X^{t). 
We also note that the definition of stopping lines can be extended to 0* by projection. 

Now, for every x e S", one can define a probability measure P* on il.* in the following way: 

— Initially, X^{0) = x. 

— The individuals on the spine move according to the strong Markov process X and die 
at the rate mi{y)R{y), when at the point y ^ S'. 

— When an individual on the spine dies at the point ?/ G ff, it leaves k offspring at the point 
where it has died, with probability {mi{x))~^ kq{x , ■) (this is also called the size-biased 
distribution of q{x, •)^). 

— Amongst those offspring, the next individual on the spine is chosen uniformly. This 
individual repeats the behaviour of its parent (started at the point y). 

— The other offspring initiate independent branching Markov processes according to the 
law P^, independently of the spine. 

This decomposition first appeared in [26]. We now have 

Lemma 3.1 (Many-to-one). Let ^ he a stopping line, such that P* -almost surely, there exists 
t ^ 0, such that {^t,t) £ Denote this time by T. Let Y he a random variable of the form 

where an -measurable random variable for every ueU . Then 

W^Ye^o R^M^))''^^M*))<^^ =w\y ^ Yu . (3.2) 

(•u,t)e^ 



{lem : many_to_one} 



Proofs of this result can be found for fixed time in [49], [36] or [38]. The proofs in [36] and 
[38] can be extended to stopping lines once the martingales that appear in the proof are still 
uniformly integrable when stopped at the stopping line Adapting the arguments of [48] or 
[11] to the continuous-time setting, one sees that this is true by the hypothesis we have placed 
on This hypothesis is also referred to as the stopping line ^ being dissecting. 

Often, we will use a simpler version of the Many-to-one lemma, which is the following 

{lem : many_to_oiie_siir 

*If R{x) is bounded from above, which will always be the case in this paper, this individual exists with 
probability one. 

^The size-biased distribution of the Dirac-mass at is again the Dirac-mass at 
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Lemma 3.2 (Simple Many-to-one). Let T = T{X) be a stopping time for the strong Markov 
process X which satisfies P^(T < cx)) = 1 for every x e S". Let / : — > [0, oo) be measurable. 
Then we have 



{u,t)e^T 



The next lemma tells us about second moments of sums of the previous type. To state it, 
we define for a stopping time T for X, the density kernel of the branching Markov process 
before .^t, by 



Pt[x 



,dy,t) = B''^ 2 l(x„(t)edj/, t<T(x„)) 
ue,_V{t) 



(3.3) 



Lemma 3.3. Let H be the hitting time functional of a closed set F <^ S on D{£') which 
satisfies P'^{H < oo) = 1 for every xe S. Let f : ^ [0,oo) be measurable. Then we have 



{eq: def ^density} 



-Clem : many_to_two> 



( 2 /(^"(*)) 



E^ 



\ 2 {f{Xu{t))f 



r ( pH{x,dy,t)R{y)m2{y)(By[ ^ 



dt (3.4) {eq : many_to_two} 



This lemma can be proven using the Many-to-few lemma from [38] (which is valid for 
stopping lines as well by the same argument as the one above) or with Lemma 3.1, by noting 
that 



{u,t)e^H {u,t)eSeH (",t)e^H {^,s)e^H, v^u 

For an intuitive explanation of the terms appearing in (3.4), see the proof of Proposition 18 
in [10]. 

Taking for X the space-time process (li,t)j>o of ^ possibly non-homogeneous strong 
Markov process (lt)i>o with paths in D{^) and the closed set F = <§ y. {t}, for some t ^ 0, 
we obtain the following useful corollary, which appeared already in [62] and [66] in the homo- 
geneous case. 

Lemma 3.4. Let f : x [0, go) be measurable and let t ^ 0. Then we have 



{lem : raany_to_two_f ix 



2-1 



e(-,o) ^ 2 /(y4t),t))^J =e(-'0)[ 2 ifiYu{t),t)f 

+ ££p(x,dy,s)i?(y,t)m2(y,t)(E(^'^)[ ^ 



ne,/r(t) 



ds (3.5) {eq : many_to_two_f i 



3.4 Doob transforms 

As in the previous subsection, we assume for simplicity that the strong Markov process X 
admits a representation as a conservative strong Markov process X with paths in D{S), 
which is killed at a rate R{x), where R : S' ^ [0, cc) is measurable. Let H be the hitting time 



{sec : doob} 
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functional of a closed set F c on D{S'). Furthermore, \et h : F ^ [0, 1] be a measurable 
function. We can then extend the function h{x) to S by setting 

h{x) = w[ W KXuit)) 

We are going to assume that h{x) > for all x e S'\F. Then for all such x we can define a 
law on Q by 

where the multiplication is in the sense of a Radon-Nikodym derivative. Now define 
Q{x) = Y,Q{^)h{x)'-\ and fc) = 

By (3.1), we now have (dropping the symbol for better reading and setting H = H{X0)) 

k 

/i(x)P^(da;) = P"(dX)(l(H<^)/j(X(/7))g(X(F),A;)nP^(^)(do.«) 

+ li^^^H)h{X{C-)fq{X{C-),k) fl P^(^-)(du;«) KXu{t))) . 

i=l (ti,t)e^H(i.j{')) 

If we denote by X^ the process X stopped at i/, and the law of X^ under by (P^)^, 
then the last equation and the strong Markov property give 



h{x)Pl{du:) = {PT{dX'')[l^H<oKX{H)) + l^^^H)h{X{C-))Q{X{C-)) ^ 

. ^ \ (*^'^} {eqibrnp conditioned 

x(l(^<^)P^W(dc.(0'^)) + l(^^^)g,(X(C-),fe)nPf ^'^(dc^^^^))- 



In particular, integrating over k, uj^'^\ i = 1, 2, . . ., and X0{t) for t e [H, <^), we get that 

h{x) = (ET[l(^H<oHXm + l(^^^H)h{X{C-))Q{X{C-))). 
We can therefore define a law on the paths in D{(o) stopped at H by 

Pl{dX) = (/i(x))-i(l(^(.Y)<c)/i(^(^)) +l(c^//(x))/i(^(C-))Q(^(C-))) X (PTidx), 

where the multiplication is again in the sense of a Radon-Nikodym derivative. Then (3.6) 
yields the following decomposition of the law P^: 

— As long as a particle has not hit the set F yet, it moves according to the law P^^, and, 
when it gets killed at the point y, spawns k offspring according to the law qh{y, which 
initiate independent branching Markov processes according to the law P^. 

— When a particle hits the set F at the point y, it continues as a branching Markov process 
according to the law P^. 
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If R{x) = R, one gets a simpler characterization of the law P^: In this case, h{x) is a harmonic 
function for the law of the stopped process under P^, whence we can define the Doob 
transform 

P^idX) = (/i(x))-i(l(H.„) + l^H<oo)h{X{H)))p'^{dX^). 

Then the law P^ is obtained from the law Pj^ by killing the process at the time-dependent 
rate RQ{x)l(^t<H)- 



4 Branching Brownian motion with absorption at a critical hne 

From this section on, q{k) will denote a law on {0, 1, 2, . . .}. We define m = ~ 
and m2 = Xifc ^(^ ~ l)*?!^) suppose that m > and m2 < oo. We let cr > and define 
Co = a\^2m. In this section we are studying the branching Markov process where, starting 
with a single particle at the origin, particles move according to Brownian motion with variance 
0"^ > and drift —cq and branch at rate 1 into k particles according to the reproduction law 
q{k). At the point —x, we add an absorbing barrier to the process, i.e. particles hitting this 
barrier are instantly killed. Formally, we are considering the process up to the stopping line 
•^H^y ) where H-y is the hitting time functional of the point —y. It is well-known since Kesten 
[44] that this process gets extinct almost surely. As a consequence, the number of particles 
absorbed at the barrier, i.e. the random variable 

Zy = #^H-y, 

is almost surely finite. By the strong branching property and the translational invariance 
of Brownian motion, one sees that the process {Zy)y^Q is a continuous-time Galton-Watson 
process, a fact which was first noticed by Neveu [59] (see [7], Chapter III or [39], Chapter 
V for an introduction to continuous-time Galton-Watson processes). Neveu also stated that 
the infinitesimal generating function u{s) of this process has the representation u = ip' o 
where ■0 is a so-called traveling wave of the FKPP (Fisher-Kolmogorov-Petrovskii-Piskounov) 
equation: Write f{s) = Y^k Then is a solution of the equation 



{sec : critical_line} 



^" - CqiIj' = 11) - f oijj, 



(4.1) 



{eq: travel ing_ wave} 



with '(/'(— oo) = 1 and ip{+a:)) is the extinction probability of the process, i.e. the smaller root 
of f{s) = s. For a proof of these results, see [54], Section 3. 

In the same paper [59], Neveu introduced his multiplicative martingales, which he used to 
derive the Seneta-Heyde norming for the martingale e~'^"^Zj^. He proved that in the case of 
binary branching, one has 



coye-'^yZy 



(4.2) 



as y ^ CO, where > almost surely. His proof relied on a known asymptotic for the 
traveling wave ip, namely that 



l-tPi-x) ~ J^xe-^o^, 



as x 



00, 



(4.3) 



for some constant K > 0. It was recently shown [67] that this asymptotic is true if and only 
if £'[Llog^ L] < 00 and the proof of (4.2) works in this case as well. We also still have in this 
case, for every x e R, 

E[e--''^''W] = iP{x), (4.4) 



{eq:def_Wy} 



{eq: tw_asymptotics} 



{eq: W_laplace} 
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{lera:y_zeta> 



a fact which was already proven by Neveu [59] for dyadic branching. 

In [10], further properties of the hmit W have been estabhshed under the hypothesis of 
dyadic branching, namely 

P(W > X) ~ — , as X ^ OD, (4.5) {eqiW^tam 

X 

and 

E[14^1(^y^^-)] — log X ^ C4.6, as X ^ CX), (4.6) {eq:W_expec> 

for some constant 04,6 e R. Equation (4.5) has been proven in Propositions 27 and 40 of [10], 
and (4.6) appears in the proof of Proposition 39 of the same paper. Their arguments were 
very ingenious but indirect and although they could be extended to general reproduction laws 
with finite variance, we will reprove them here directly under (probably) minimal assumptions, 
based on methods of [54]. The main result in this section is 

-[prop : W} 

Proposition 4.1. //£'[Llog^L] < 00, then (4.5) holds. If E[L log^ L] < 00, then (4.6) holds. 

See also [24] for a proof of (4.5) in the case of branching random walk. Before proving this 
result in the next subsection, we state a lemma which is immediate from (4.2) and the fact 
that Zy is almost surely finite (see also Corollary 25 in [10]): 

Lemma 4.2. Suppose £^[Llog^L] < 00. For any ri> 0, there exist y and such that 

y^T]-^ and P{\Wy - W\ > r]) + P{^h^^ ^ U x ^ V- 

4.1 Proof of Proposition 4.1 

Rescaling space by Cq"^ , we can always come down to the case cq = 1, which we will assume 
to hold for the rest of the proof. Define x{^) = -E[e~^^]. Our first result is: 

Lemma 4.3. Suppose that -E[Llog^L] < go. Then, 

X"W - \, (4.7) 
as X ^ 0. If furthermore £'[Llog'^ L] < 00, then 

X"(A) = '"(-^)l (4-8) ■(eq:chi_asymptotic_2 

A 

where r(A) ^ and r(A) dA < 00. 

Proof. Define </>(x) = 1 — Tp{—x), such that u{s) = (l)'{(l)^^{s)). By (4.1), we have 

^Ct)"{x)+^'{X)= f{l-^{x))-{l-cj){x)). (4.9) <e,:pM> 

Then by (4.3), we have 

(j){x) ~ Kxe~^, as X ^ CC. (4-10) <;eq:phi.asymptotic} 

Setting g{s) = 2[/(l - s) - 1 + f'{l)s] ^ and p = + 0', we get from (4.9), 

p{x) = —p{x)+g{(p{x)). (4.11) {eq:rho} 



■Clem: chi} 



{eq: chi_asymptotic> 
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As in the proof of Theorem 1.1 in [54], we will study the function p through the integral 
equation corresponding to (4.11), namely 

= e-^' U(0) + e^(7(</'(2/)) = e"^' p(0) + TT^ . (4.12) {e,:„egrai> 

Now, by Theorem B of [13] (see also Theorem 8.1.8 in [14]) we have for every d ^ 0, 
loe'^ - log'' - 1 J 

— ^c/(s)ds<oo <^ — ^g(s)ds<oo ^[L log^+"' L] < oo. (4.13) {eq:bingham} 

Jo Jo 

Furthermore, by Proposition 3.2 in [54], we have —u{s) ~ s, as s — > 0, and by (4.10), we have 
e?^ ~ (log l/s)/s, as s ^ 0. Under the hypothesis £'[Llog^ L] < cc, we therefore have 

ds < 00, 

Jo -Hs) 

whence, by (4.12), 

p{x) ~ Ke"^, as X ^ 00, (4-14) {eq:rho,asymptotic} 

where the constant K is actually the same as the one in (4.10), see the proof of Theorem 1.1 
in [54]. Now, from (4.4), we get xW = 1 - 0(-logA), whence, by (4.11) and (4.12), 



"(A) = -^p'(-logA) = - + ^(-Aj ^ ^eyg{cl,{y))dy -g{cP{- log X)) 

r eV(y)5'(0(y))dy), 



(4.15) {eq:chi_rho> 



A A 

where the last equation follows from integration by parts. This proves (4.7), with the constant 
K instead of 1, since the last integral vanishes as A — > 0. Now, setting 



(^) = -t( r eV(y)<7'('A(y))dy 

^ ^ J- log A 



'-log 

we first remark that r(A) ^ 0, since the integrand is negative for y e R. By the Fubini-Tonelli 
theorem, we then have 

rr(A)dA= C-\ r eV(y)5'(<A(y))dydA 

Jo Jo J- log A 



I 

Jo 





e^-'''yU-Hy)9'{y)dy, 

Jo 



which is finite if and only if S[Llog^L] < oo, by (4.13) and the fact that e*^ ^^^''(f>^^{y) ~ 
(log^ l/s)/s. This proves (4.8), again with the constant K instead of 1. 

The previous arguments worked for every traveling wave ■0. In order to show that that 
the constant K is equal to 1 in our case, we use Neveu's multiplicative martingale (see also 
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[51], Theorem 2.5). We know [59, 25], ((1 — (/){x + y))^^)y^o is a martingale for every x G R, 
bounded by 1. By (4.2) and (4.10), we get by dominated convergence, for every x e R, 

x(Ke^) = lim E[e-^y^''''^y] = lim ^[(1 - (Piy - x)fy] = 1 - (/)(-x) = xie""). 

y—i-co y— >co 
This yields K = 1. □ 
Remark 4.4. Choosing arbitrary initial points xo,xi e R instead of in (4.12), one sees that 

J -00 rco 
eyg{^{y)) dy = e^V(xi) + e^5(0(y)) dy. 

In particular, since p is bounded, letting xq — > — oo and xi +oo yields 



■CO 

eyg{^{y))dy = l. 



One could hope (see the proof of Proposition 4.1 below) that this helps in determining the 
constant C4,6, but apparently this does not seem to be the case. 

Proof of Proposition ^.1. We define the function 

Vn{x) = rx"P(P^ e dx) = E[W\w^x)l 
Jo 

such that with x^") denoting the n-th derivative of x, we have for A > 0, 

x(")(A) = (-!)" e-^MK(x). 
Jo 

If £'[Llog^L] < 00, Proposition 4.1 and Karamata's Tauberian theorem ([33], Theorem 
XIII. 5. 2 or [14], Theorem 1.7.1) now yields 

V2(x) ~ X, as X ^ 00. (4.16) {eq:V2} 

By an integration by parts argument (see also [33], Theorem VIII. 9. 2 or [14], Theorem 8.1.2), 
we get (4.5). Now suppose that i?[Llog^L] < co. By Lemma 4.3, we have x'{^) ~ ^og \ 
c G R, as A ^ 0. By Theorem 3.9.1 from [14] (with ^{x) = 1), this yields 

Vi(x) — logx ^ 7 — c, as X ^ 00, 

where 7 is the Euler-Mascheroni constant. This is exactly (4.6). 

□ 

5 Branching Brownian motion in an interval 

{sec : interval} 

In this section we study branching Brownian motion killed upon exiting an interval. Most ideas 
in this section (except for Section 5.4) stem from Sections 2 and 3 of [10] and for completeness, 
we will reprove some of their results with streamlined proofs. However, we will also extend 
their results to the case of Brownian motion with variable drift. 
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5.1 Notation 

During the rest of the paper, the symbol C stands for a positive constant, which may only 
depend on the reproduction law except in Section 7, where it may also depend on some 
other constants which will be specified. Its value may change from line to line. If a subscript 
is present, then this subscript is the number of the equation where this constant appears for 
the first time (example: 05,23). In this case, this constant is fixed after its value has been 
chosen in the corresponding equation. If X is any mathematical expression, then the symbol 
0{X) stands for a term whose absolute value is bounded by C|X|. 

Recall the definition of q{k)^ m and m2 from Section 4 and the hypotheses on m and 
m2- From this section on, we further define cq = -slim. Furthermore, in this section, we let 
a ^ vr/co and set 



2m 



From (5.1), one easily gets the basic estimate 



^ Co - /i ^ 



TT 

2/ia2 



(5.1) 



(5.2) 



We then denote by the law of the branching Markov process where, starting with a 
single particle at the point a; e R, particles move according to Brownian motion with variance 
1 and drift — /x and branch at rate 1 into k particles according to the reproduction law q{k\ 
Expectation with respect to is denoted by E^. On the space of continuous functions from 
R+ to R, we define and Ha to be the hitting time functionals of and a. We further set 
H = Hq a Ha- Then note that the density kernel of the branching Brownian motion below 
^H, as defined in (3.3), has a density with respect to Lebesgue measure given for t > and 
x,ye (0,a) by 



where was defined in (2.4). 

Now, let / e '^^(R^o, R^so) be non-decreasing, with /(O) 
called a barrier function. We set 



(5.3) 



0. Such a function will be 



max 



Now define 



IJit = H + 



1 



(5.4) 



(5.5) 



such that /iQ = ^ and //f ^ /i for all t ^ 0. We denote by P^ the law of the branching Brownian 
motion described above, but with infinitesimal drift —fit- Expectation with respect to P^; is 

denoted by Ej and the density of the process is denoted by p{{x,y). 

The above definitions can be extended to arbitrary initial configurations of particles dis- 
tributed according to a counting measure z/ on (0, a). In this case the superscript x is replaced 
by v or simply omitted if is known from the context. 



{sec : interval_notati 



{eqimu} 



{eq: cO_mu} 



{eq: density_bbm} 



{eq : def _f „norm} 



-teq:mu_t> 



5.2 The processes Zt and Yt 

Recall from Section 3 that the set of particles alive at time t is denoted by ^{t). We define 

^,a(t) = {u e J^{t) : H{Xu) > t}, 
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where H was defined in the previous subsection. Now set w{x) = ae^^^ sin(7r2;/a) and 
define 

Zt= 2 ^i^uit)) and Yt= 2 e'^(^"W-'^). 

Then Zt is a martingale under P^, since e^^w{Bt) is a martingale for a Brownian motion with 
drift — /i killed at and a, which is easily seen by Ito's formula, for example. Furthermore, it 
is easy to see as well that Zt is a supermartingale under P^. 

The following lemma relates the density of BBM with variable drift to BBM with fixed 
drift: 



-Clem:mu_t_density} 



Lemma 5.1. 



where I Err I < 



p{(x,y)=pt(x,y)e-<=«/(*/'^')+E-, 



1 _1_ * _1_ TT' 



Proof. By the Many-to-one lemma and Girsanov's theorem, we have 



exp ( - J ds^ W^^ (exp ^ fM^ - l-^dB,) , Bt e dy, H > t ] . 



By integration by parts, we have 



fig dBs = fitBt - fJ^Bo - Bs di^ig. 
Jo Jo 



(5.6) 



(5.7) 



-teq:005} 



{eq: ipp_mu_B> 



Since Bt e (0, a) for all t ^ 0, we have 









Bsdfis 




Jo 


Jo 





d^^iiriioo-. 



(51 



{eq:010} 



Furthermore, 



such that 



Finally, 



nt/a 



2\2 



2 a 



f 

Jo 



2 



t - fif{t/a' 



2a4 ' 

2a4 
ll/'lloo 



HtBt-uBt = ^nt/a^]B, ^ 
Equations (5.3), (5.6), (5.7), (5.8), (5.9) and (5.10) now give 

p/(x,y) = pi(x,y)e-^^(*/'^')+'^- 
and the lemma now follows from (5.2). 



(5.9) {eq:020} 

(5.10) 



-[eq:030> 



□ 
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(5.12) 



{eq: Zt^variance} 



{prop : quantities} 

Proposition 5.2. Under any initial configuration of particles, for every t ^ 0, we have 

Bf[Zt] = Zoe-""-^(*/'^')+^", (5.11) {e,:Zt_expectation> 

and if in addition jj, ^ cq/2, then 

Ye.vj{Zt)^Ce'^"(^^Zo + Yo). 
Furthermore, we have for every t ^ ( without hypothesis on ), 

and for t ^ , 

E/[yt] ^ Ce^"^. (5.14) M:Yt> 

Here, Err has the same meaning as in Lemma 5.1. 



(5.13) {eqiYt.weal!} 



Proof. Equation (5.11) follows from Lemma 5.1 and the fact that Zt is a martingale under 
P^. In order to show (5.13) and (5.14), it suffices by Lemma 5.1 to consider the case without 
variable drift. We first suppose that t ^ a?. By (5.3) and (2.8), we get 

W[Yt\ ^ e^(^-") r e^^pt{x, y) dy ^ Ce^(^-«) sin(7rx/a) - sin(7ry/a) dy. 
Jo Jo a 

The last integral is independent of a. Summing over x yields (5.14) as well as (5.13) in the 
case t ^ . Now, if t < a^, by the Many-to-one lemma and Girsanov's theorem, we have 

W[Yt] = e*VFf^[e^(^*~'^\ HQAHa<t]= e^'*/(2a2)^x^^^ a Ha < t]e'^(^^'*). 

Summing over x yields (5.13). 

Throughout the proof of (5.12), we use a constant C, which depends only on the repro- 
duction law q{k) and which may change from line to line. By Lemma 3.4, 



[Z2] = E^[ 2 w{Xu{t)f\+2m2[ f p{(3;,2/)(Ej(._^^)[Zt_,])2dsdy. (5.15) M:bb._2nctao.ent> 



By Lemma 5.1 and the fact that Zt is a martingale with respect to the law P^, (5.15) yields 



E^ 



■[Z?]^Ce^^"M 2 w{X^{t)f]+j''j\,{x,y)w{y)^dsdy 



(5.16) -Ceq:050> 



Now we have for x e (0, a), 

w{xf = (asin(7r2;/a)e-^("-^'))2 ^ TT^{a - xfe'^^"^"-^^ ^ Ce^(^~"), 
because /i ^ cq/2 by hypothesis. This yields 



Si:=E-[ 2 MXu{t)f 
ueJ^o,a{t) 



(5.17) 



{eq:Sl} 
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by (5.13). Now, by (5.3) and (2.13), we have 

32-.= r f Ps{x,y)w{y)Usdy = ae^'^''-''^ C ae'^'^y-'^') sin^Try/ a) r{x,y,t) dy. 
Jo Jo Jo 

Lemma 2.1 now gives 

^2 ^ Cae^^""'"-^ J e^^y s'm^{TTy/a)(tsm{TTx/a) sin(7ry/a) + ayj dy 

° . rCO (5-18) {eq:S2> 

^ Cae'^^^-") ( sm(7rx/a) ^ + - ) e'^'yy^ dy, 
\ a-^ aJ Jo 

the last hne following again from the change of variables y ^ a — y and the inequality sin x ^ x. 
Using again the fact that /i ^ co/2, equations (5.16), (5.17) and (5.18) now imply 

E^[Z|] ^ Ce^^" + e^'^'^-'^A. (5.19) ^e,:Zt.2nd.on,ent> 



If we write the positions of the initial particles as xi, . . . then by the independence of 
their contributions to Zt and by (5.19), 

Var;(Zi) = 2 ^ 2e^:[^'] ^ e3E"2E-'[z2] ^ Ce^^- (Lz^ + , (5.20) 
by (5.19). This proves (5.12). □ 
5.3 The particles hitting the right border 

In this section we recall some formulas from [10] about the number of particles hitting the 
right border of the interval. We reprove these formulae here for completeness and because 
Lemma 2.1 makes their proofs straightforward. For most formulae we will assume that / = 0, 
i.e. that we are working under the measure P. Only Lemma 5.6 contains an upper bound on 
the expected number of particles for general /, which will be useful in Section 7. 

For a measurable subset 5 i= R, define Rs to be the number of particles killed at the right 
border during the (time) interval 5, i.e. 

Rs = #{(n,t) : u e ^{t) and Hq{Xu) > Ha{Xu) = teS}. 

The following lemma gives exact formulae of the expectation and the second moment of Rs- 

{lem : R_moments } 

Lemma 5.3. For every x e (0,a), we have 

E^[i?5] = e^(^-")r(x,5), (5.21) {e,:RS.x_expe 

na rco 2 

E^[ii|] = E^[i?5] +m2e'^(^""M dye^'^y-'^U dtef^^p^{x,y)I''{y,S - (5.22) 

Jo Jo 

We will first prove a more general result, which will be needed in Section 6.4. 

{lera : R_raany_to_one} 

Lemma 5.4. For every x e (0,a) and any measurable function f : R_|- R_|-, we have 



c} 



-Ceq : RS_x_2ndinoment> 



[ 2 /(t)l(x4i)^a)] =e'^(--'^)£/(^)/'^ 



E^ 



(x, ds). 
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Proof. Recall that Hq and Ha denote the hitting time functionals of and a and H = Hq a Ha- 
Then note that W!^^{H < ao) = I for all x e [0, a]. We then have 



/(i)l(X„(t)=a) 



e^^^f{Ha)l^Ho>H. 



by Lemma 3.2 



= e^'^''-''^WS\e^"''fiHa)l^Ho>Ha^t)] by Girsanov's transform 

= eM(^-a) r*j(s)/a(2.^ds) by (2.12). 

Jo 

□ 

Proof of Lemma 5.3. Equation (5.21) follows from Lemma 5.3 and (2.14) by taking / = I5. 
Equation (5.22) follows from Lemma 3.3 and (5.21). □ 

Lemma 5.5. For any initial configuration u, and any ^ s ^ t, we have 

irlt — s) , ^ / „ , , s , •^s 



s,t] 



Zo\ ^ C5.23 ( Yo A E,/a^{l A (t - s)/a^)Zo ) , 



(5.23) {eq : Rt_nu_expec} 



where Eg is defined in (2.7). Furthermore, if fi ^ Co/2 and ^ t ^ a^, then for each x e (0, a), 

t 



WRl ^ ^5,24 I -^w{x^ 



_|_ ^K^-a) 



(5.24) 



-[eq:Rt.nu_variEince} 



Proof. We have Ei?^ = ^u(dx)E^Rt, such that (5.23) follows from (5.21) and Lemma 2.1. For 
the second moment, we have by (5.22), 

E^[i?2 - Rt] = mae^^^-'^) f dy e'^'^y"'^ f ds ef^^p^^ix, y)P{y, t - s)^ 

Jo Jo 

^m2e^(^-") f dye^'^y-''^I''{y,tfr{x,y,t) 
Jo 

^ Ce'^^^-") r dy e^(^-") {t/a^ sin(7ry/o) + 1)^ 
Jo 

X (at/a^ sin(7r2;/a) sin(7ry/a) +a^^{x a y){a — {x v ?/))) 

Performing the change of variables y 1-^ a — y in the integral and making use of the inequalities 
a~^{x A y){a — {x V y)) ^ a — y and sinx ^ x, we get 

rOD 

E^[i?2 _ R^] ^ 0(l)e^(^-") ism{TTx/a)t/a^ + 1) dy e^^^^iy + yH/a^ + yH^/a^) 

Jo 

^ Ce'^(^"")(sin(7r2;/a)t/a2 + 1)(1 + t'^/a^), 

for some constant C, which does not depend on /i by the hypothesis /i ^ co/2. Using the 
hypothesis t ^ and (5.23) yields (5.24). □ 

Lemma 5.6. Let f be a function as in Section 5.1. Then for every x e (0,a), we have 

WjiRs] ^ W[Rs\. 
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Proof. As in the proof of Lemma 5.3, we have 



Ha ,,2, 



by Girsanov's theorem. Now, we have by (5.7), on the event [Hq > Ha}, 

rHa rHa 

jjLt dBt = iJ.{a -x) + a{nHa - A*) - Bt dfit > fi{a - x), 
Jo Jo 



{Ho>HaeS) 



since Bt ^ [0, a] for t e [0, Ha]- This gives 
by the proof of Lemma 5.3. 



^{Ho>HaeS) 



□ 



We finish this section with a lemma which hnks BBM with absorption at a critical line to 
our BBM with selection model. 

Lemma 5.7. Let ^ 1, y ^ 1, ^ cq/2 and f be a barrier function (defined in Section 5.1). 
Suppose that s/a ^ y + C ^"^^^ ll/ll ^ V^- -^^^ collection of space-time points 

with 

Xi = a-y + {cq- fi)ti - f{s/a^), z = 1, . . . ,iV, 
and ^ C for all i. Define Z = Y,i w{xi), Y = Yji e^"^^'-"^ and Wy = coye~'"'yN. Then, 

In particular, for large a, we have 

Y ^ Z/y. 

Proof. By (5.2) and the hypotheses /i ^ cq/2 and C ^ 1, we have for all i, 

fCn - 

Xi = a-y + 0( 

Hence, by (5.2) and the hypotheses /i ^ co/2, ||/|| ^ ^/a and ^/a ^ y + Ci 

e^X. _ g/.a-COS/ (l^0(^^y (5.25) {eq:250> 

Furthermore, since x — x^/3 ^ sinx ^ x for x ^ 0, and by the hypotheses y ^ 1, and a ^ y + C 

sin(7rxj/a) = sin(7r(a — Xi)/a) = — y ^1 + O^— . (5.26) {eq:270} 

The lemma now follows by summing over (5.25) and (5.26). □ 



{lem: critical_line} 
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5.4 Penalizing the particles hitting the right border 

In this section, let {Uu)ueU be iid random variables, uniformly distributed on (0, 1), indepen- 
dent of the branching Brownian motion. Furthermore, let p{t) : R+ e (0, 1] be measurable 
and such that p{t) = for large enough t. Recall that H = Hq a Ha- We define the event 

E = {t{u,t) e ■■ Xu{t) = a and Uu ^ p{t)}. 

Our goal in this section is to describe the law = P^(-|i?). We first note that 



{sec :weakly_ conditio 



P^(da;,£;)=P^(da.) W (l(x4i)^a)+p(i)l(x„(*)-a) 

(u,t)eLt 

In order to apply the results from Section 3.4, we define 

h{x,t) = P(^'*)(E) 

CO 



fc=0 



qix,t,k) = q{k)h{x,tf''^ /Q{x,t) 



(5.27) {eq:180> 

(5.28) {eq : def_ gamma} 

(5.29) {eq:def3eta> 

(5.30) -teq:def_btilde} 



By the results from Section 3.4, under the law P^', the BBM stopped at is the branching 
Markov process where 

— particles move according to the Doob transform of Brownian motion with drift — /x, 
stopped at and a, by the space-time harmonic function h{x,t), and 

— a particle located at the point x at time t branches at rate Q{x,t)lxe{o,a): throwing k 
offspring with probability q{x,t,k) and 

— a particle located at or a does not branch. 

We have the following useful Many-to-one lemma for the conditioned process stopped at 
the stopping line = -^H ^ i- Define the function 



e(x,t) =Y,k{l- h(x,tf-^)q{k) ^ 7712(1 - h{x,t)). 



(5.31) {eq:def.e} 



Lemma 5.8. For any measurable function f : [0, a] — > R+, we have 



■Clem : many_to_one_coii 



e4 2 fiXuit)) 



/i(x,0) 



(5.32) {eq : many_to_one_cond 



In particular, if we denote by p{x,y,t) the density of the P^-BBM, then 

Proof. By Lemma 3.2 and the description of the law P^ given above, we have 

E-[ 2 f{Xu{t))\ =t^^^[/(X^^,i)^(^//Ai,i^ At)eSo''"-(^-^)Q(^-^)d^ 

where m{x,t) = Xifc(^' ~ ^)Qix,t,k). Applying Girsanov's transform yields (5.32). Equation 
(5.33) follows from (5.32) applied to the Dirac Delta-function f = 6y, y e (0, a), together with 
(5.3). □ 



(5.33) {eq:ptilde_estimate} 
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The previous lemma immediately gives an upper bound for the quantities we are interested 



m: 



Corollary 5.9. For any x e (0, a) and t ^ 0, we have 

B^'lZt] {h{x,0)y^'E''[Zt], (5.34) {eq:conditioned_Z_up 

£"^[^1] ^ (/1(X,0))"^E^'[Z|], (5.35) {eq:conditioned_Z2} 

E^'lYt] {h{x,0)y^E'=[Yt]. (5.36) {eq:conditioned_Y> 

Proof. Equations (5.34) and (5.36) immediately follow from (5.33) and the fact that 7(2/, t) ^ 1 
for all y and t. In order to prove (5.35), we note that by Lemma 3.4 and the description of 
the conditioned process, 

E-[z2] = e4 2 wiXuit))^] + f r pix,y,t)m2iy,t)Qiy,t) (b^-^'^IZ])' dydt, 

ueN{t) J Jo Jo ^ ^ 

where fn2{x,t) = Xlfc>o^(^ ~ ^)Q{x,t, k). Equation (5.35) now follows from (5.15) together 
with (5.33), (5.30) and (5.34). □ 

The following lemma gives a good lower bound on E^[Zt]. We define 

pt = sup p{s). 

se[0,t] 

-[lem: conditioned_Z_l 

Lemma 5.10. Suppose /i ^ co/2 and ^ i ^ a'^ and p{s) = for all s ^ t. We have 

Ef[Zt\ > W{X) (l - PtC5,37{t/a^ + (1 A (a - X)"^)) . (5.37) {eq:conditioned_Z_lo 



In order to prove it, we will need the following estimate on h{x, 0), which will be sharpened 
in Lemma 6.2. 

Lemma 5.11. 

1 - h{x, 0) ^ pti7rw{x)t/a^ + Ce'^^^"")) ^ ptCe''^'^^''^ (^(a - x)t/a^ + l) . 
Proof. By Markov's inequality, we have 

1 - h{x,0) = P"(#{(n, s)e^t: Xu{s) = a, Uu ^ p{t)} ^ 1) 
^ B'^miu, s)e^t- X^{s) = a, C/„ ^ p{t)}) 

The lemma now follows from Lemma 5.5 and the inequality sinx ^ x, x e [0, vr]. □ 

Proof of Lemma 5.10. Since p{s) = for all s ^ t by hypothesis, we have h{y,t) = 1 for all 
y e (0,a). Lemma 5.8 and the second property of the Brownian taboo process (see Section 
2.3) now imply 

B-[Zt]^w{x)W,l^,,[e-^o<^^'')<i^'^^uj{x)[l-W,^^^^^^ (5.38) <eq:..o> 



-[lem:h_estimate} 
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by the inequality e ^ ^ 1 — x for x ^ 0. By (5.31), Lemma 5.11 and the hypotheses, we have 
for every y e (0, a) and s ^ 0, 

e{a -y,s)^ e{a - y, 0) ^ ptCe'^'y^yt/a^ + 1) ^ PtCe-(^«/3)^. 

By Lemma 2.2 and the fact that the law of the Brownian taboo process is preserved under 
the map y i— > a — y, this gives 

rt 



^'^taboo 



^ e{Xs,s)ds s£:PtC{t/a^ + (iii{a-x)^. 



The lemma now follows from (5.38) and (5.39). 

Finally, we study the law of Rt under the new probability. 
Lemma 5.12. We have for every x e [0, a], 

E-[Rt] -ptB'^iRh ^ E-[i?t] ^ iKx,0))-^E-[Rtl 
and if there is a p e [0, 1], such that p{s) = p for s ^ t, then we even have 

E"[i2i] ^ E[Rt]. 

Proof. Let be the stopping line 

m = {{u,s)e^Ha --s^t}. 
We have by definition of the law P, 



(5.39) {eq:220} 
□ 



{lem: conditioned_Rt} 



(5.40) {eq: conditioned_Rt} 
^5.41) {eq: conditioned_Rt_c 



E-[Rt] 



Now the denominator is h{x,0) by (5.28), which yields the right-hand side of (5.40). The 
left-hand side follows by noticing that 

B-[Rt Yl {l-p{Xu{s)))]>E-[Rt{l-ptf^]^-E-[Rt]-ptB-[R^t]. 
For (5.41), we note that if p{s) = p for s ^ t, then by (5.42), 



(5.42) {eq:230> 



Since (1 — p)^ is decreasing in /c, this yields (5.41). 



□ 



6 The system before a breakout 

In this section, we are studying branching Brownian motion with drift — /i and absorption at 
until a breakout occurs, an event which will be defined in Section 6.1 and which corresponds 
to a particle going far to the right and spawning a big number of descendants. In (6.20), 
we decompose the system into a particle conditioned to break out at a specific time T (this 
particle will be called the fugitive) and the remaining particles, which are conditioned not to 
break out before time T. These two parts will be studied seperately, the former in Section 
6.4 and the latter in Section 6.3. Before that, in Section 6.2, we study the law of the time 
of the first breakout, showing that it is approximately exponentially distributed. First of all, 
however, we start with the necessary definitions: 



-[sec :bef ore_breakout 
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6.1 Definitions 

We will introduce several parameters which will be used during the rest of the paper. The 
two most important parameters are a and A, which are both large positive constants. The 
meaning of a is as in the previous sections: It is the right border of an interval in which the 
particles are staying most of the time, and a breakout will be defined below as the event that 
a particle hits a and then spawns many descendants. The parameter A has a more subtle 
meaning and controls the number of particles of the system and with it the intensity at which 
particles hit the point a. In Section 7, we will indeed choose the initial conditions such that 
Zq % Ke^, where k is a fixed constant. 

In [10], the parameter a was called (which we changed for typographical reasons) and 
a and A were related by the equation 



{sec :bef ore_breakout 



a = — f log + 3 log log iV - A 
Co V 

where was a parameter representing the approximate number of individuals in the system. 
The parameter A then represented a shift of the right barrier. Although this choice of param- 
eters may be more intuitive then ours, we found it technically more convenient to drop the 
parameter A^ altogether, and work only with a and A instead. 

As in [10], when we study the system when a and A are large, we will first let a go to 
infinity, then A. Thus, the statement "For large A and a we have. . ." means: "There exists 
Aq and a function uq^A), such that for A ^ Aq and a ^ ao(^) we have. . .". Likewise, the 
statement "As A and a go to infinity. . ." means "For all A there exists ao{A) such that as A 
goes to infinity and a ^ ao{A). . .". We further introduce the notation o^(l), which stands for 
a deterministic term independent of the initial conditions of the process and which goes to 
as A and a go to infinity. Furthermore, o(l) will denote a term which goes to as a goes to 
infinity (with A fixed). 

The remaining parameters we introduce are all going to depend on A, but not on a. First of 
all, there is the small parameter e, which controls the intensity of the breakouts. Indeed, when 
Zq ^ Ke"^, the mean time one has to wait for a breakout will be approximately proportional 
in e. Morally, one could choose e such that e"^/^ « e « A^^, but for technical reasons we 
will require that 



e ^ Ce.i^"-^^, and 



(6.1) {eq:ep,upper> 

^6.2) {eq : ep_lower} 



Another protagonist is ry, which we will choose as small as we need and which will be used to 
bound the probability of very improbable events, as well as the contribution of the variable 
Y. It will be enough to require that 

-2^, (6.3) 



r] ^ e 



which, by (6.2), implies 



12 



(6.4) 



The last parameters are y and C, which are defined as in Lemma 4.2, with rj there being the 
?7 defined above. Note that the parameters r], y and appeared already in [10] and had the 
same meaning there. 



{eq: eta} 



{eq:eta_ep} 
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We can now proceed to the definition of tlie process. Recall the definition of /i in (5.1). 
We will always suppose that a is large enough, such that 



/i ^ Co/2. 



(6.5) 



-teq:mu_cO} 



As in Section 5.1, we denote by P^' and the law and expectation of branching Brownian 
motion with drift — starting from a particle at the point j; e R, and extend this defintion 
to general initial distributions of particles according to a counting measure u. Recall from 
Section 3.1 that ^{t) denotes the set of individuals alive at time t. We want to absorb the 
particles at and do this formally by setting 

yKo{t) = {ue^{t):Ho{Xu)>t}, 

where Hq is again the hitting time functional of 0. 

Instead of absorbing particles at a, we are now going to classify them into tiers in the 
following way: Particles that have never hit the point a form the particles of tier 0. As soon 
as a particle hits a it advances to tier 1. Say this happens at a time tq. In order to advance 
to tier 2, a particle has to come back to the critical line a — y + (cq — /x)(t — tq) and then back 
to a again. Here, y is a large positive constant to be defined later. 

Formally, let u e U, t ^ 0. We define two sequences of random times (T„(n))„^_i and 
((Tn(it))n>0 by r_i(n) = 0, cro{u) = and for n ^ 0: 



Tn{u) = inf{s ^ (7n{u) : Xu{s) = a}, 
<7„+i(ti) = inf{s ^ Tn{u) : Xu(s) = a - y + {cq - fi){s - r„(ti))}, 

where we set inf = oo. We now define for i ^ the stopping lines 

^f-* = {(n, s) eU x[0,t]: s = ti{u) and ne J^q{s)}, I ^ -1, and 
= {(u,s) e [/ X [0,t] : s = a/(n) and ne ^(s)}, / ^ 0. 

That means, contains the particles of tier / at the moment at which they touch the right 
barrier and ^/'^ contains the particles of tier / at the moment at which they come back to the 



^6.6^ {eq: def _tau_sigma> 



(6.7) 



critical line. Note that the sets 
for every 1^0. We also set 



Pj^''' and ^/'^ are increasing in t and 



{eq: def _Rline} 



(6.8) -(eq:def_Lline} 



In order to extend the definitions of the variables Zt and Yt to the current setup, we could 
simply replace .AQ,a{t) by .A(){t) in their definition (see Section 5.1). However, it will be more 
useful to take special care of the individuals u for which ti{u) ^ t < (Tj+i('u) for some I ^ 0, 
since these are in some kind of "intermediary" state which is difficult to analyse. We therefore 
define the stopping lines 

^yf^^'-* = {(u, s) e U X R+ : u e ^(s), r/_i(ii) ^ s < r;(u), 

and either t < (Ji{u) = s or (Ti{u) ^ t = s}, / ^ 0, (6.9) {eq:def_Niine_i} 



and 



(6.10) {eq:def_Nline} 



1^0 
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such that c/^^^'' contains the particles of tier I that have already come back to the critical line 
at time t, as well as the descendants of those that haven't, at the moment at which they hit 
the critical line. We then define for I ^ (recall that w{x) = asm{7Tx/a)e^^^~°'^), 

For any symbol S and ^ /c ^ /, we define, 



5(fc;0 = ^ = 5'a;oo)^ g ^ ^(o+)^ 



i=k 

For a particle (n, s) G .^^^'^ , we define the stopping line 

_ |(^^^^^ e f/ X R+ : v e ^(r), < (7;,r) and r = o-i+i(t;)}. 

This stopping line yields a collection (X„(r),r — s)^^^ ^^^^y(u,a) of space-time points, and we 
denote by y("'^) and W^"'"^ the corresponding quantities from Lemma 5.7. Of course, 

we have chosen the stopping line in such a way that the variable Wy^'^^ follows the same law 
as the variable Wy defined in (4.2). We also define r^x^ = max^^ (r — s). We then 

define the event 

^(M = |^(M>^eAj^|^(n,s)><^|, (6.11) u„ 

which is called the event of a breakout, since e measures the number of descendants of 

the particle (u, s) (the inclusion of the "bad" event {t^-ax > C} is for technical reasons). The 
particle u is then also called the fugitive. We set 

PB = P"(5(0'°)), (6.12) <e^-„y 

and define the law of BBM started at a with the first particle conditioned not to break out: 



Q«(.) = P'^(. |B^) 



1 -PB 



where we set B = b'^^'^I We further set Z = and Wy = W^^'^^ and note that by 

Lemmas 5.7 and 4.2, we have for large a, 

P-(|Z - ^W\ > 2n) + V%T^^^ ^0<V, (6.13) 
Co 

where W is defined as in (4.2). Hence, by (6.3) and (4.5), we get 



{eq:Z.U} 



PB = (— +Oa{1) +0{1))^, (6.14) {eq:pB> 

V Co / ee^ 



which goes to as ^ and a go to infinity, by (6.2). Furthermore, (6.13) yields for large A and 



Q'^iZ] = (E[-m(^H/^,,A(i+„(i))+o(^))] +0(r?ee^))(l + 0(pb)) 

^ (6.15) {eq:QaZ> 

= —{A + loge + C6.15 + oa{1) + o(l)), 
Co 
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by (4.6), (6.1), (6.3) and (6.14). In particular, we have for 1 and large a, 

Q°[Z] ^ CA. (6.16) 
Moreover, by (4.5), (6.13) and (6.3), we have for A^ 1 and large a, 

Qa|-^2j ^ Cee^. (6.17) <;eq:QaZsquared} 

We now define for every / e Z^o the time of the first breakout of a particle of tier I, 

tW(l^) = inf{t ^ : w e (J S^"'")}, (6.18) {eq:def^Ti> 



and set 



min tW. (6.19) 



We denote by the fugitive of the breakout that happened at time T. 

Now fix t > 0. We want to describe the system conditioned on T = t. For this, suppose 
that at time the particles are distributed according to a counting measure v = Xi^^i'^^^i- 
Define pi = F'^{i < \ T = t), which yields a law (pi)^^^ on the initial particles, depending on 
1^ and t. Since the variable T, the time of the first breakout, is the minimum of the variables 
Tj, i = 1, . . . , n, the times of the first breakout of the BBM descending from the particle i, we 
can decompose the process into 

n n 

( Yl d^^'^ \t = ^) =YiPi (P'^Hdw^'^ \T = t)xYl P""' {dJ^'^ \T>t)y (6.20) <eq: decomposition} 

That is, we first choose according to the law {pi)'^^-^ the initial particle that is going to cause 
the breakout. This particle spawns a BBM conditioned to break out at time t. The remaining 
particles spawn independent BBM conditioned not to break out before time t. 



6.2 The time of the first brealcout 

We want to prove that the random variable T defined above is approximately exponentially 
distributed with parameter p^vrZo/a^, which is the statement of the following proposition: 

Proposition 6.1. Let ^ t ^ a^/(3C6.38^) o^nd suppose that Yq ^ e^^. Define 9 = pbttZq. 
Then, for A and a large enough, we have 

3/1 , /i+/„3 



■Csec : time_breakout} 



-[prop:T> 



P{T >t) = exp\^-et/a''\^l + 0{At/a-' +pB)j +0{pBYo)j. (6.21) {eq:T> 

The proof proceeds by a sequence of lemmas. Lemma 6.2 gives a estimate on P(t('^) > t). 
This is used in Lemma 6.4, in order to obtain an estimate on P"(r > t), using a recursive 
argument. Finally, Proposition 6.1 is proven by combining Lemmas 6.2 and 6.4. 

-ClemiTi 

Lemma 6.2. Let ^ t ^ a^. Define 9 as in Proposition 6.1. Suppose that pB ^ 1/2 and 
Yq ^ e^^. Then, 

P(r(0) > i) = exp - 9t/a^{l + 0{pb)) + ©(pb^^o))- (6.22) <eq:To; 



36 



Before proving Lemma 6.2, we establish a weaker estimate on P[T^^^ > t)). 

Lemma 6.3. Let ^ t ^ a^. Define 6 as in Proposition 6.1. Suppose that ps ^ 1/2 and 
^ 1- Then 



-Clem:TO_ small} 



P(r(0) > t) = exp ( - et/a3(l + 0{pb{1 + ^ot/a'))) + 0(pb>o))) • 
Proof. We have for t ^ 0, 



E 



^1 -PB 



(0) 



(6.23) -[eq:TO_small> 



(6.24) {eq:30S} 



since by the strong branching property, the random variables Z^"'^) are independent condi- 
tioned on By Lemma 5.5, and the assumption t ^ a^, we have 



and 



|E[i?f ] - TrZot/a^\ ^ C5.23Y0, and 



E[(i?fV] ^ 2{7rZot/a^f + C{TrZot/a^ + Yq), (6.26) {eq:3iB> 

where the last line follows from the inequality {x + y)'^ ^ 2(x^ + y^) and the assumption 
Iq ^ 1. By Jensen's inequality and (6.25), we have 



(6.25) {eq:310> 



E 



^1 -PB 



(0) 



^ E 



,log(l-pB)E[i?f'] 



exp ( - et/a^ + O {pB{Ot/a^ + Yq)) ) , (6.27) {eq:3m 



since | log(l — x)\ ^ x + for x ^ 1/2. This gives the lower bound in (6.23). For the upper 
bound, equations (6.25) and (6.26) together with the inequality (1— p)" ^ l — np+n{n—l)p'^/2 
give 



E 



(1 - pb)^' ' ^ 1 - et/a'' + pBC{et/a:' + Fo) + {Ot/a:')'\ 



(0) 



3^2 



(6.28) {eq:325> 



The lemma now follows from (6.24), (6.27) and (6.28) together with the inequality 1 — x ^ 
e-^. □ 

Proof of Lemma 6.2. Let xi, . . . ,x„ be the positions of the initial particles. Since the initial 
particles spawn independent branching Brownian motions, we have 



p(tW >t) = ]^p^>(t(°) > t). 



(6.29) {eq:340} 



Define Zi = w{xi) and yi = e^^^* '^K Then trivially yj ^ Yq for all i, and therefore, since 
lo ^ e by assumption. 



/i(a - Xj) = I logyil ^ I loglol > 1, for all i. 
As a consequence, by the inequality sinx ^ x for x ^ 0, we have 

Zi = ayism.{'KXi/a) ^ -Kfi^^e^^^"'^''^^^ n{a — Xj) ^ -Kfi^^Yo\ loglol. 



(6.30) •teq:345> 



since the function x xe ^ is decreasing for x ^ 1. By Lemma 6.3, (6.5), (6.29) and (6.30) 
and the hypothesis t ^ a^, we now have, 

P(t(0) >t) = Yle,,p(^-pB7Tz,t/a^(l + 0{pB{l + Yo\logYo\))') +0{pBy^)). 

i 

Since Yq \ logYol ^ 1 by hypothesis, this proves the lemma. □ 
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In the following lemma, note that according to the definition of the tiers, a particle starting 
at a starts immediately in tier 1. 



Lemma 6.4. Let ^ t ^ a^/(3C6.38^)- Then, for large A and a, 

Q'^(T >t)^exp(^- pBTT^Q^iZ] (1 + 0{A^ + pb)) + 0(r?) 
Proof. We have 

Q«(r>t) = Q"[ Yl P^"^'^(r>t-s)] ^Q°[P''(r>t)], 

where u = s)ey^^^ ^Xu{s)- Since T > t implies T^''^ > t, we have 

p^(r >t) = p^iT > t|r(°) > t)p'^(r(°) > t). 



(6.31) 



(6.32) 



{lem:T.a} 



•(eq:T_a> 



-Ceq:360} 



(6.33) 



-Ceq:362} 



Let Z and Y be as in Lemma 5.7. By the definition of Q" and Lemma 5.7, we have Q"-almost 
surely Y ^ rjee"^, such that for large A, y ^ e^^ and pbY ^ 7]C by (6.3) and (6.14). By 
Lemma 6.2, we now have for large A, 

P^(t(°) > t) ^ exp - 9Z (1 + 0(pb)) + 0(r/)) . (6.34) ^eq:363> 

As for the first factor in (6.33), we have, with the notation from Section 5.4, with p{s) = pB, 



p'^(r > t|r(°) > t) = p'^(r > t) = p'^J Yl Q"(r>t-s) 

By Jensen's inequality and (5.41), this implies 

P'^iT > t|T(0) >t)^ Q°(T > t)E''[«r] ^ Q-(r > t)E-[/?f ] 
Now, by (5.23), we have, by Lemma 5.7 and y^^ ^ rj, 



(0) 



ri 



(6.35) 



(6.36) 



-Ceq:364} 



{eq:365> 



Equations (6.33), (6.34), (6.35) and (6.36), together with Jensen's inequality, now imply 

Q"[P'^(r > t)] ^ Q"(r > t)Q"m(-Va^+^) X exp ( - eQ'^[Z] (1 + 0{pb)) + 0(r/)) . (6.37) <eq:37o> 
Now, by (6.16), (6.32) and (6.37), we have, 

Q''{T>t)^-^ ^exp( -9Q''[Z]{1 + 0{pb)) + 0{i])), (6.38) ie^.sr.y 



with 5 = C6.38(^i/a^ + 7/(1 + j4)). By (6.3) and the hypothesis on t, we have 6 ^ 1/2 for large 
A, whence (1 — 6)^^ ^1 + 26. Raising both sides in (6.38) to the power (1 — 5)^^ yields the 
lemma. □ 
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E[i?f = vrt/a^Zo + 0{Yo 



Proof of Proposition 6.1. We have the trivial upper bound P(T > t) ^ p(T(o) > t), and 
Lemma 6.2 now imphes the upper bound in (6.21). For the lower bound, we note that as in 
the proof of Lemma 6.4, we have by Jensen's inequality and (5.41), 

p(r > t) = p(r > 1 1 r(°) > t)p(r(°) >t)^ ct{T > t)E[i?^°*]p(2-(o) > ^y 39) ^^^^3^^^ 

By Lemma 6.4, and since At/a^ = 0(1) by hypothesis, we have 

Q''(r >t)^ exp (0{pBAt/a^ + rj) 

and by (5.23), we have 

such that, since At/a^ = 0(1) by hypothesis, 

Q"(r > t)E[«f''] ^ exp (O {9t/a^{At/a^ + ti/pb) +v + YoiPB + v)) ) (6.40) <e,:3B2> 

The lower bound in (6.21) now follows from (6.39), (6.40) and (6.22), together with the fact 
that 7] ^ Cp% by (6.3) and (6.14). □ 

{lem : moment s_T} 

Lemma 6.5. Define 9 as in Proposition 6.1. Suppose that Yq ^ e^"^ and let Q and n e N. 

Then, for large A, 

n I 

E[(r/a3 + a)"l(r^,,,)] ^: E[(r(o)/a3 + a)"l(^(o)^,3)] ^ C ^ J^^^^y^ (6-41) <e,:x._> 

Furthermore, if ^ 13 = oa{A^^), then for large A and a, 

E[(r/a3) 1(^/^3)^^] = 9-\l + 0{AP + pb)) + 0((/3 + 0-i)e-^(^^)) (6.42) <e,:T^e.pec> 

Proof. We first note that we have, for n ^ and 7 > 0, 



fc=0 

Now, we have 



(t + a)"e-^*dt= V ■ . (6.43) 



•(eq:398} 



E[(T(0)/a3 + a)"l(^(o)^„3)] =['(* + a)" P(T(o)/a3 e dt) 

^ n r(t + a)"-i P(t(°) > ta^) dt + a". 
Jo 

The second inequality of (6.41) now follows from Lemma 6.2 and (6.43), since pB ^ as 
A goes to infinity. The first inequality follows in the same way, using the trivial fact that 
P(r > to^) ^ P(r(0) >ta?). For the second part, we note that 

E[(r/a3)l(r/,3)^^] = P(r > ta^) dt - /3P(r > f3a^), 
Jo 

and by Proposition 6.1 and the hypothesis on Yq, we have for t ^ P and large A and a, 

P(r> to^) = {i + 0{pB))exp{-9t{l + 0{A/3+pB))). 
Equation (6.42) now follows from the last two equations. □ 



39 



We now show how we can couple the variable T with an exponentially distributed variable: 

{lem: coupling.!} 

Lemma 6.6. Suppose there exists a universal constant k, such that c^^Zq = k + 0{e^^'^) and 
that Yq ^ rjZQ. Then there exists a coupling (T,V), such that T is a (V) -measurable and the 
random variable V is exponentially distributed with parameter pbg'^ttk and such that for large 
A and large a, we have P(i?coupl) ^ Ce'^, where i?coupl is the event 

Scoupi = {\T/a^ -V\> £3/2} u {|T/(aV) - 1| > V^) ^ e^}. 

Proof. For brevity, we define 7 := pbc^ttk. Let F be the tail distribution function of T, i.e. 
F{t) := P(T ^ t). It is clear that T has no atoms except 00. We can therefore define a random 
variable U which is uniformly distributed on (0, 1) by setting 

U = F(r)l(r<„) + U'F{cc)l^T=a,), 

where U' is a uniformly distributed random variable on (0, 1), independent of T. Now we 
define V = —7"^ log U. Then V is exponentially distributed with parameter 7 and T = 
F^^{e^'^^), where F^^ denotes the generalized right-continuous inverse of F. Hence, T is 
(T(y)-measurable. On {T < 00}, we have by Proposition 6.1, for a large enough, 

V = -j-HpBe''7re-^ZoT/a\l + 0{AT/a'^ +Pb))) + 0{pbYo) ^ ^ 

„ „,„ „ , (6.44) {eq:390} 

= T/a^{l + 0{e^/^ + AT/a^ +Pb)) + 0{KpBe'^r]), 
by the hypotheses on Zq and Yq. Hence, by (6.2), (6.4) and (6.14), we have for a large enough, 

\T/a^ -V\= 0{e^/^T/a^ + A{T/a^f) + 0{e^l^). 
But now we have by Lemma 6.2, for large A and a, 

P(r/a3 > e^I^/VA) ^ P(r(0)/a3 > £3/4/^) ^ ^^-o(e-VVVZ) ^ ^2/2^ (6.45) 
by (6.2). Furthermore, we get from (6.44), 

V/{T/a^) = 1 + 0(6^/2 ^ AT/a^) + {0{eri) + o{l)) / {T / a^) , (6.46) m:394> 

and by (6.4), we have by Proposition 6.1, 

P(T/a3 ^ Ver?) = 0(??/Ve) ^ £^2, (6.47) {eq:396> 

for large A and a. Equations (6.44), (6.45), (6.46) and (6.47) now prove the lemma. □ 

6.3 The particles that do not participate in the breakout 

{sec :bef ore_breakout 

In this section, we fix t ^ /{?>CQ,'i^A). We are going to study the system conditioned not to 
break out until time i, the law and expectation of which are denoted as in Section 5.4 by P 
and E, respectively, hence 

p(-) = p(-|r > t). 

Under the law P, the process stopped at a t then follows the law P from Section 5.4, with 

P(S) = PB'^is^t) + (1 - PBWiT ^t-s), (6.48) {eq:Cef_ps> 
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such that by Lemma 6.4, (6.3) and (6.16), for large A and a, 

Pt ^ CpB- 



(6.49) 



{eqipbar^estimate} 



As in the proof of Lemma 5.12, one can then show that Q"[.Z^] = (1 + 0{pb))Q,"'[Z] and 
Q''[Z2] ^ {1 + 0{pb))Q''[Z^], such that by (6.15), (6.16) and (6.17), we have for large A and 
a, 



Q"[Z] = —{A + loge + C6 15 + oa{1)), Q"[^] ^ CA, and Q"[Z^] ^ Cee^ (6.50) {eq:Qhatz> 
Co 

We define two filtrations {'^i)i^o {-^1)1^0 t>y 
such that and cz Mi cz for every I. Now define 



(.u,s) ^(0 



(n,s)ej^f' 



and recall from Section 6.1 the definition of .^0''^ -^0^^ ^0 the corresponding 

quantities for Y. 

Lemma 6.7. Suppose Yq ^ 77^0- have for all I ^ 1, and /an/e ^ and a, 
E[zg I ^ (^ + C6.5lPB)Q"m(iT + C6.51??)^g"'^ 

/n particular, 

E[zg] ^ ((^ + C6.5lPB)Q'^m(iT + C6.5ir?))'Zo. 
/n t/ie case I = 1, we also have for large A and a, 

B[ZI^^] ^{tt- Ce.53PB)mZ]i^ - Ce.53V)Zo- 

Proof. We have 



-Clem : Zl_exp_upperbou 



Qe\z\R 



■4 ! 



^6.51) {eq : Zl_exp_upperbouii 

(6.52) {eq : Zl_exp_upperboun 

(6.53) {eq : Zl_exp_lowerboun 

(6.54) {eqiZl.Hl} 



since conditioned on M„(i-\\ ,, the random variables Z^"''^) are iid under P of the same law as 

XI At' 

Z under and independent of =^^('-1)^^) t>y the strong branching property. Now, we have 



E[i?f I 



]= 2 E^"(^)[i?£)jr>t-s]. 



(6.55) {eq:400} 



By (5.40) and (5.23), the right-hand side of (6.55) is less than or equal to 



t 



max /i(Xu(s),s) ^ 
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and we have I0 ^ Z0 ^■'/y ^ Z0 ^■'r/, P-almost surely, by the definition of the event F 
in Lemma 5.7 for / ^ 2 and by hypothesis for / = 1. Furthermore, by Lemma 5.11, (6.5) 
and (6.49), we have (/i(x,0))~^ ^ 1 + CpB, as soon as ps is small enough. Combining these 
inequalities with (6.54) and (6.55) gives (6.51). Equation (6.52) follows easily from (6.51). 
Now, in the case / = 1, we have % = ^0 by definition. Let u = Xil^i denote the initial 
configuration. By (5.40) and (5.23) and (5.24), we have 



?(o) 



(0)1 



?(o)^2^ 



^ TT^Zo - C5.23>0 - PtQ 



5.241, 



This yields (6.53), since Yq ^ r]ZQ by hypothesis. 



□ 



In applications of Lemma 6.7, we will often sum the right-hand side of (6.52) over all / ^ 
We therefore define tg.se to be the solution of 

(vr + C6.5iPB)Q°[Z](i6.56/a^ + Ce.5iv) = 1/2- 



(6.56) -[eq:def_tcrit} 



We now turn to the variance of Z0 in the cases that are of interest to us, namely, for 
/ = 1,2. 



Lemma 6.8. We have for I ^ 1 and large A and a, 



-(lem: Zl_variance} 



VarfZ, 



(0 



^,)^Cee^zl^-'\t/a^ 



v)- 



In particular, we have for 1=1, 



Var(zJJ^) ^ Cee^'Zoit/a^ + rj). 



Proof. We have 



Var(zg 



1) = E[Var(zg | ^) | + Var(E[Z^^ | J^] | ^^,_i) 

= VarQ„(Z)E[i2f-') | + q,^[Z]Y^T{B^-^^ \ % 



(0 



(6.57) -Ceq:430> 



By Lemma 5.5, the assumption t ^ c? and the fact that Yq ^ Z^^jy on Gq (in the case / = 1), 
we have 



E[i2f~^^ I ^ C(t/a^ + r?)Z^ "^ and 

The lemma now follows from these equations, together with (6.2) and (6.50). 
Lemma 6.9. Suppose that Zq ^ Ce^ . Then we have for large A and a, 

^^(zg^) ^ Ce^Zo(e^(t/a3 + r^f + £^2(^/^3 ^ ^^3 ^ ^4(^/^3 ^ ^^4^^ 
Proof. We have by repeated application of Lemmas 6.7 and 6.8, 
V£rr(zg)) = Vrr(E[zg) | ^1]) + E[Va:r(zg) | ^^i)] 



Y^i{Rt'^ I ^^,_i) ^ C{t/a' + v)zt'l 



□ 



{lera : Z2_variance> 



^ E[{CA{t/a^ + ?7)ZJ!V] + E[Cee'^Z^S>{t/a^ 



^ CA\t/a^ + 7?)2(Va:r(zg)) + E[Z^^^f) + Cee^{t/a^ + 77)E[Z«] 
^ C{A'^ee'^{t/a^ + r]fZo + A^{t/a^ + 7])^Z^ + ee^A{t/a^ + v?Zo). 
The hypothesis on Zq now proves the lemma. 



□ 
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■(lem:Zl_t} 

Lemma 6.10. We have for all I ^ and for large A and a, 

P(|zi') - > K I ^0 ^ CK~^Z^^ + l(4)>c/./p,)- 
In particular, suppose that t ^ ie.SG; then 

P(|zf - >K)^ CK-^Zo + Ce-^Zo/{Ke). 

Proof. We decompose 



P(|zf - 4^1 > 2Ani^0 ^ P(|Z(') - E[Z« |^;]| > I^^O + l(|E[z«|^#,]-z(')|>/r)- (6.58) « 



Now, we have by the conditional Chebychev inequahty: 



P(|zi') - E[zf I m > K I ^ (6.59) M:.e.> 



By (5.35) and (5.19), 



Yai{Zf^\^l)^ 2 E(^"^'^''^[^f^] ^ C'Zg. (6.60) ^e,:463> 

As for the second term in (6.58), we have by (5.11), (5.34), (6.49) and Lemmas 5.10 and 5.11, 

\E[ZP I ^l] - I ^ CpbZ^§ , (6.61) M:4S5> 

since xe~^^ ^ C by (6.5). Equations (6.58), (6.60) and (6.61) now finish the proof of the first 
inequahty. The second inequality follows readily by taking expectations and using Lemma 6.7 
and Markov's inequality. □ 

{lem:Z3.and.Y} 

Lemma 6.11. Suppose that t ^ ie.se- Then for large A and a, 

P(Zf+^ > A') ^ CK"^A^(t/a^ +7?)^Zo, and (6.62) <eq:Z3_estimate} 

P{Yt > AT) ^ CK-^r]Zo. (6.63) {eq:Y_estimate} 

Proof. First note that we have /i(x, 0) ^ 1/2 for large A. We now have by (5.34) and (5.11), 

E[ZP] ^ 2E[zg)] ^ Cr?Zo(^ + Ce.sir?)'. 

Using the hypothesis, summing over / ^ 3 and applying Markov's inequality yields (6.62). For 
(6.63), we note that by (5.36) and (5.13) and /i(x, 0) ^ 1/2 

E[y/')] ^ CE[4')] ^ CA-SE[zg] ^ Cr^Zo2-\ 

by Lemma 6.7 and the hypothesis. Summing over / ^ and using Markov's inequality finishes 
the proof. □ 

43 



6.4 The fugitive and its family 

We now describe the BBM starting from a single particle and conditioned to break out at a 
fixed time t. We could describe this system by similar methods as those employed in Section 
3.4, but since we are only interested in first moment estimates, it is faster to use the Many- 
to-one lemma instead, which is the method of the proof of the following lemma: 

Lemma 6.12. Let f : [0, a] — > R+ be measurable, ^ r ^ a'^ and t ^ t. Let p : R_|_ [0, 1] 
be measurable with p{s) = for s ^ t and CQ^ePt ^ 1/4. Denote by P the law associated to 
p(s) as in Section 5.4- Then, 



{sec : fugitive} 



{lem : raany_to_one_f ug 



?(0) 



n f/ X {r} 7^ 



^ C6.64W^taboo 



(r - s)) 



ds 



(6.64) {eq : many_to_one_f ugi 



Proof. The left-hand side in (6.64) equals 



LHS 









Ef 




l(Ha(X„)edr) 





(6.65) 



By Lemma 3.1, the numerator of the right-hand side of (6.65) equals 

According to the description of the conditioned process in Section 5.4 and the description of the 
spine in Section 3.3, the particles on the spine spawn on average Q{x, s) k{k — l)q{x, s, k) ds 
particles during an interval [s, s + ds], which is less than or equal to m2h{x, s) ds. Conditioning 
on the trajectory of the spine and using (5.33) now yields 



{eqiSOO} 



NUM ^ m2e'^''W' 



,(r - s)) 



dsl 



{Ho>HaedT) 



Applying Girsanov's theorem to this expression and Lemma 5.8 to the denumerator in (6.65), 
we get 



LHS ^ 



m2W, 



killed 



Vo^''iI.M.M^JiX.ir-s))]ds 



where VFj^^^^lJ is the law of a bridge from x to a of length r of a Brownian motion killed 
at and a. But since the taboo process is obtained from the killed BM by a space-time 
Doob transform, this is the same as W^^^'^^. As in the proof of Lemma 5.10, we have, by the 
hypotheses on r, 

K^oll^-^^^''^''-'^'''] > 1 - Ce.em- (6.66) 



This implies the lemma, by (6.49) and the hypothesis on pf. 



□ 
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We now set up the important definitions. Recall that denotes the fugitive. Define 
^{i) = {u e ^o(i) : K^) a ('^.T) e f/ x (J[az(^), t«(^))}, 

^{t) = {{u,s) e (J^«(t) : a (^,r) e [/ x tK^)))} 

and 

jY{t) = {u e ^(t) : A (^,T) 6 \J x J [Ti_i(^), a/C^))}. 
We then define 

ue.jV{t) ueJ^{t) 

and 

ue.yf{t) ueJv{t) 

Note that on the event T = T*^^\ we have ^{T) = by definition. For the other particles, 
we have: 

■[lem: cousins_TO} 

Lemma 6.13. Suppose that t ^ tg.se '^"-'^ C^^Q^pt ^ 1/4. Then, 

W[Zt\T = T^^^ =t\^CA, W[Yt\T = T^^^ =t\^C, W[Rt\T = T^^'^ = t\ ^ C 

Proof. We have for every s ^ 0, E^[zi'^^] = w{x) ^ 7r(a — x)e^^^'^''^^ by (5.11). Furthermore, 
E^[yi°^] ^ Cw{x)/a ^ Ce-^'("-^') by (5.13). Finally, we have by Lemma 5.5, W[Rf^] ^ 
C((a — x) + l)e^^("~^') for s ^ a? . By Lemmas 6.12 and 2.2, we now have 

E^[zf ) I T = r(o) = t] ^ c, E^ [y/°^ I T = r(°) = ^ c, E^[i?f ^ | r = r(°) = t] ^ c. 

From the estimate on Rf\ it follows that W[Z^^^ \ T = T^^) = i] ^ CQ°[Z] ^ CA, by (6.50). 

Hence, by Lemma 6.7 and the hypothesis, we have E^[Z0^^ | T = T^*^) = t] ^ CA. By (5.34) 
and (5.11), we now have 

E^[zf+^|T = T(°) = t] ^CA, 
and by (5.36) and (5.13), we have 

E^[ya+) I ^ ^(0) _ ^ 

Since = Zj'^^ + zl:^^\ Yt = Y^'^'^ + Y^^^^ and 7/ ^ A^"^, this implies the lemma. □ 
On the event T = T^^^ , the situation is more complex, as shown by the following lemma. 

-(lem: cousins_Tl} 

Lemma 6.14. Suppose that t ^ tg.se ^^c^ CQ,QQPt ^ 1/4. Then, 

W[Zt I T = r(^) = ^ c^, w\Yt I T = r(^) = t] ^ c, E^[i2t I T = r(^) =t\^C. 

Moreover, on the event T = T^^^\ we have Z^^ ^ ee^, and 

E^[zf I T = r(i) =t\^ Cee^A{t/a^ + r?), 
E^[yj I T = r(^) = ^ Cer]e^, 
P^p(i) - > K I r = r(i) = ^ Cee^K-\t/a^ + rj) + C{t/a^ + r])/K. 
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Proof. On the event T = T^^\ conditioning on ri(^), we get two independent pieces of the 
process, one starting at x conditioned to hit a at Ti(^), the second starting at a, following 
the law of Q"^ conditioned on T = T^^^ = t — Ti(^) =: t' . Stopping this process at the line 
one of those particles then spawns BBM conditioned on T = T^^^ = t' and the others 
spawn BBM conditioned on T > t'. Now, Zt and Yt are the sums of the respective variables 
corresponding to the two pieces and the inequalities on their expectations now follows from 
Lemma 6.13. 

On the event T = T^^\ for / ^ 1, we can generalize the above decomposition and condi- 
tioning on ri('^), . . . ,Ti_i(^) we get / independent pieces of the process. On this event, we 
note that ^ ee^, since no breakout occurred before the time t = ti{^). This immediately 

gives the estimates on the first and second moment of ^0^. 

For the proof of the remaining inequalities, we note that we have by Lemma 6.7, 

00 

1=1 

by the hypothesis on t. The last three equations now follow from these results as in Lemmas 
6.10 and 6.11. □ 

Define TC'™) = min^^j^^ and r('+) = miui^^T^. 
Lemma 6.15. Suppose that Cie^ ^ Zq ^ C2e^ and Yq ^ tjZq. Then for large A, 

p(r(i+) < rW) ^ CeA, 

and 

P(r(2+) <T^^'^))^C{eAf. 

Proof. Let to := ^6.56 a /{ZC^.^fiA). By (6.2), (6.14) and (6.50), we have for large A, 

P(r(°) > to) ^ exp(-CA/e) ^ e. (6.67) {eq:777} 

Now, for the rest of the proof, let t ^ to- We have by the decomposition (6.20) of the process 
conditioned on T = t, 

n 



{lem:Tl_T2} 



pi.(^(i+) > ^ I ^(0) = t) = 2 piP^-^-i (t(^+) > 1 1 r(°) > t) p^'(r(i+) > 1 1 r(°) = t) 

i=l 



by Jensen's inequality. By Lemma 5.5 and the hypothesis on Yq, we have E'^ "^"^i [i?^'^''] ^ 
^"[R?^] ^ C{t/a^ + T])Zo. By Lemma 6.13, we have for large A, E^'[i?f^ | r(0) = t] ^ C. In 
total, we get by Lemma 6.4, for t ^ to, 



pv^^(i+) > i| j(0) _ ^) ^ (^_c(^ApB{t/a^ + rifZo + t/a^{'nZo + ApB) + v + ri^Zo)y (6.68) {e,:62o> 

By (6.67), (6.68), and the inequality 1 — ^ x, we get 

p-(r(i+) < r(0)) ^ C(E[(T/a3 + ^)H^t^,^,s)]ApbZo 

+ E[T/a=^l(r^toa3)](r/Zo + Apb) + v + rfZo + e). (6.69) 
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The result now follows by Lemma 6.5, together with the hypothesis on Zq, (6.2), (6.3) and 

(6.14) . For the second part of the lemma, we first note that by (6.20) and the union bound, 

n 

pu (^(2+) ^ ^ I ^(1) > j.io) =t) = Y, Pi (p'-^^' (t(2+) < 1 1 r(O'i) > t) 

+ p^' (r(2+) < 1 1 t(i) > r(°) = t)y 

Now we have 

pu~s,^ (r(2+) < 1 1 r(°'i) > t) ^ p'^(t(2+) < 1 1 t(°'1) > t) 

by Markov's inequality. As in the proof of Lemma 6.7, we can show that for Z ^ 1, 

Qa^^a + l) I > i I J^^] ^ C7^'(t/a3 + 7?)'i?f \ (6.70) ^e,:631> 

since we have, as in (6.50), for every / ^ 0, Q"[Z | r(i;') > t] = (1 + 0(pb))Q"[^] ^ CA, by 

(6.15) . With (5.23) and (6.49), this gives 

(^(2+) ^ ^ I ^(0;i) ^ ^) ^ C76.7iPiJ^^(i/a^ + r/)^^o. (6.71) U^-.em 

Moreover, we have 

px^y(2+) ^ ^ I ^(1) ^ y(o) _ ^) ^ ptE^[i?f ^ I r(i) > t(°) = t] 

< PBCQ''[Zf{t/a^ + v?'E''[Rf'> \ T^o) = t] by (6.70), (6.49) 

^ PBCA^{t/a^ + rjf by Lemma 6.13. 

In total, this gives 

pu^r^{2+) ^ ^ I ^(1) ^ ^(0) ^ ^) ^ CpB{A\t/a^ + rjfZo + A'^{t/a^ + T]f) (6.72) {e,:66o> 

Moreover, we have 

n 

p.(j.(2+) ^ ^ I ^(0) ^ ^(1) ^ ^) ^ 2 |^p.-5., (^(2+) ^ ^ I ^(0;1) > ^) 

+ P^>(T(2 + ) < t I T(°) > r(^) = t)) . (6.73) {e,:661> 

The first term in (6.73) has been bounded in (6.71). For the second term, we note that we 
have 

px^^(2+) ^ ^ I ^(0) > ^(1) _ ^) ^ ptC7(Q'^[i?f ^ I r(^) > t]E^[i?J°^ I r(°) > r(^) = t] 

and by Lemma 6.14, together with (6.2), (6.14), (6.49) and the hypotheses on Zq and Yq, we 
get 

P^(r(2+) < 1 1 r(°) > r(i) =t)^ Cps (^A\t/a^ + r,f + ee^A{t/a^ + r?)^) ^ CA{t/a^ + 

(6.74) {eq:663} 
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Equations (6.71), (6.73) and (6.74) now yield 

pu^j.{2+) ^ ^ I j,{0) > ^ ^) ^ CpBA\t/a^ + 7]fZo + CA{t/a^ + r^f, (6.75) i.^:.m 
and (6.72) and (6.75) then yield 

pi. (^(2+) ^ ^ I ^(0;i) _ ^) ^ CpBA^{t/a^ + 7]fZo + CA{t/a^ + r]f. (6.76) {e,:668> 

The second part of the lemma now follows from (6.76), by integrating over t from to to ^-nd 
using Lemma 6.5 and (6.67). □ 

7 The system with the moving barrier 

{sec :moving_barrier} 

We will now define properly the BBM with the moving barrier. We will still use all the defini- 
tions from Section 6.1, with one notational change: Recall that by (6.20), we can decompose 
the process into two parts: the first part consisting of the particles spawned by the ancestor 
of the fugitive, and the second part consisting of the remaining particles. As in Section 6.4, 
the quantities which refer to the particles of the first part will be denoted by a bar (e.g. Z) or 
check (e.g. Z). The quantities of the second part will be denoted with a hat in this section 
(e.g. Z), in reference to the law P from Section 6.3. Furthermore assume from now on that 
there is a constant k, such that for each A and a large enough the initial distribution satisfies 
\e~^ZQ — k;| ^ e'^/^ and Yq ^ rfZ^. The constant k will be regarded as universal, in the sense 
that the terms denoted by 0(), oaO and o(l) may depend on k. 

Suppose further that we are given a family {fx)x^o of non-decreasing functions fx G 
^^(R, R+), such that for each x ^ 0, fx{t) = for t ^ 0, /a;(+co) = x and for each 
6 > there exist Mx = Mx{5), Mt = Mt{5), such that 

• Mx{S) CO as 6 ^ 0, 

• ll/xll ^ for all X e [0,M^], and 

• fxit) ^ X - 6 for all t ^ Mt, 

where ||/|| is defined in (5.4). It is easy to construct such a family: Take any non-decreasing 
function / e 'Tf^(Ti+,'R+), such that f{t) = for all t ^ and f{t) = 1 for all t ^ 1 and 
define fx = xf for x ^ 0. Then \\fx\\ ^ H/IK^^ v x^), whence this family satisfies the above 
conditions with Mx{S) = {\\f\\S)-^/^ a {\\f\\d)-^ and Mt{6) = 1. 

Now suppose we are given a BBM with constant drift —fi starting from the initial config- 
uration vq. We are going to define for each n e N define a stopping time T„ and a barrier 
process (Xi^"^)te[T„_i,T„] as follows: 

1. We set To = and xj^^ = 0. 

2. Denote by T the time of the first breakout of the BBM absorbed at and by the 
fugitive, as in Section 6.1. We set X^^^ = for t e [0, T]. 

3. Define 

1 f Zt \ 

A = log J V 1 , (7.1) {eq:Delta.def} 

Co V Ke^ ) 

where is defined in Section 6.1. 
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4. Define T[ = T + T^iP and Ti = (T + a^'l'^) v T[. Note that T[ and therefore also Ti is 
a stopping time for the BBM. Now define 

xr^ = /A((t-rO/a2),te[r,ri]. 

We then give the particles an additional drift —{d/dt)x[^^ for t e [T{, Ti], in the meaning 
of Section 5.1. 

5. We have now defined Ti and X^^^. We further define vi to be the measure formed by 
the particles at time Ti, which have never hit 0. To define T2 and , we repeat the 
above steps with the process formed by the BBM started from those particles, with the 
definitions changed such that the barrier process starts at x!^^ = time starts at 
Ti etc. 

6. We now construct the barrier process x['^^ from the pieces by x['^^ = x[^\ if t e 

Remark 7.1. The random line formed by the particles at time Ti which have never crossed the 
barrier X^^^ is not a stopping line in the sense that we have defined it, but in Jagers' sense 
(see Section 3.2), such that the strong branching property applies here as well. It is even a 
simple stopping line in the terminology of Biggins and Kyprianou [11]. 

Recall the definition of the phrase "As A and a go to infinity" from Section 6.1. Our main 
theorem is the following: 

{th: barrier} 

Theorem 7.2. As A and a go to infinity, the process {Xt)t^o = i-^to'^c'^/-K^ ~ converges 
in the sense of finite-dimensional distributions to the Levy process {Lt)t^o with Lq = and 
cumulant K^{X) given by 

rOD 

Ki^{X) = log E[e''^^^] = f A(l0g K + C) + CQ e*^^' - 1 - iAj;1(^^1) A(dx), (7.2) {eq:cuinulant,levy} 

Jo 

where A(da;) is the image of the measure x^^dx by the map x Cg ^ log(l + x) and c e R zs 
a constant depending only on the reproduction law q{k). 

A stronger convergence than convergence in the sense of finite-dimensional distributions 
is convergence in law with respect to Skorokhod's topology (see [32], Chapter 3). Obviously, 
the convergence in Theorem 7.2 does not hold in this stronger sense, because the barrier 
is continuous but the Levy process is not, and the set of continuous functions is closed in 
Skorokhod's topology. However, if we create artificial jumps, we can rectify this: 

-[th:barrier2> 

Theorem 7.3. Define Jt = X^\ ifte [T„, Tn+i), for n e N. Then as A and a go to infinity, 
the process {X^)t^Q = {Jta'^c^J-K'^ ~ ^t)t^Q converges in law with respect to Skorokhod's topology 
to the Levy process defined in the statement of Theorem 7.2. 

Define the sequence {Gn)n^-i of "good events" by G_i = Q and Gn to be the intersection 
of Gn-i with the following events: 

• i^n has support in (0, a), 

• (^U X {Tn} and r„ > (for n > 0). 



49 



• |e-^ZT„ - k| ^ £3/2 ^ ?/Zt„. 

The core of the proof of Theorems 7.2 and 7.3 will be the following proposition: 

{prop:piece> 

Proposition 7.4. Fix A e R. Suppose that P(Go) = 1. Define 70 = ti/ [c^pBe^) ■ Then there 
exists (5 > 0, such that for n ^ e^^^"^/^ and large A and a, we have P(G„) ^ 1 — ne^^^ and 

l0gE[e*^^^"] = nK"Vo(i^K(A) + iXA + 0^(1) + 0{e^)), (7.3) {eq:piece.fourier} 

with Kf^{X) defined as in Theorem 7.2 and where Oy!i(l) and 0{e^) may depend on A. 

Remark 7.5. The process Zt approximately measures the population size, in the sense that 
the number of particles at the time t + 6a^, where 5 is a small constant, is approximately 
27rc^^e'^°"o~'^ x Zt [10]. That is why we stated in the introduction that the number of particles 
stays almost constant in our model. We could indeed show for example that the process of the 
number of particles converges to a constant in the sense of finite-dimensional distributions. 
Since we plan to show much stronger results in a second work, we will not prove this here, 
however. 

7.1 Proof of Proposition 7.4 

In this subsection, we are under the hypotheses of Proposition 7.4, i.e. we suppose that uq has 
support in (0, a), \e^^ZQ — k\ ^ e"^/^ and Yq ^ r]Zo. 

The particles on the stopping line ,At- In a first step, we will describe the state of the 
system at the stopping line .y)^, defined in (6.10). Recall that this stopping line consists of 
those particles, for which ai ^ T < ti for some /, and of the descendants of those for which 
Ti ^ T < c7;-|_i for some /, as soon as they hit the critical line. This latter case applies in 
particular to the fugitive , for which T = Ti{'^) for some I. We will show in this paragraph 
that the following events happen with high probability and give first and second moment 
estimates of the quantities appearing in the definitions: 

= {\Zt- - {Z'^-''^ + Zg>)\ ^ e'/^A]} n {Zt - Z^^^'^^ ^ e^/e} n {Yr ^ e-'/'}, 

First of all, since the process of the descendants of the fugitive after time T follows the 
law P'^(- 1 [r^^'^^)^), we have for A and a large enough, 

P(GFug) ^ PsH^^iZ^^'^^ > e^/e) + P'^ir^^J^ > ()) ^ Ce\ (7.4) <e,:.4s> 

by (6.2) and (6.3). 

As for Gbuik) let A be large enough, such that y/e ^ te.se/o^^- Recall the decomposition of 
the BBM conditioned on T = t given by (6.20) and denote by xq the position of the particle 
that is the ancestor of the fugitive. We have as in the proof of Lemma 6.2, for large A, 

w{xo) ^ 7rfi~^Yo\logYo\ ^ Cr]\logr]\e'^ ^ e^/^e^, 

by (6.3), whence \e~^ZQ — k\ ^ Ic'l'^ and Yq ^ '^Z{). We then have by Lemma 6.7, for 
t ^ sfeo?, 

E[e-^zg^2) I T = = (vr + 0(pB))e-^ZoQ'^[Z](i, + 0(7?)) + O^A^^j^ + r?)^). (7.5) .e.no 
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Furthermore, by Lemmas 6.8 and 6.9 and the inequahty (x + y)"^ ^ 2(x^ + y^), we have 
Var(e-^zg^2) IT = t) ^ C{e{j,+0{r^)) + eA{j,+0{r^)f + eA\j, + 0{i^)f +A\j,+0{r^)f) 

(7.6) {eq:712> 

Lemma 6.5 and (7.5) now give for large 

[(■() {eq:714} 

= co7o(A + loge + c + OA(l)), 

by (6.1), (6.4) and (6.50). Note that by (6.14), 70 = e(l + 0^(1) + o(l)). Similarly, (7.5) and 
(7.6) and Lemma 6.5 give 

E((e-^zg'Vl(T^Via3)] = 0(^2e2^. (7.8) .e,:..s> 

Likewise, if we define 

then we have by (7.7), Lemmas 6.5, 6.10 and 6.11, the hypotheses on Zq and Yq and the union 
bound, 

P(G^ T ^ Vea^) ^ C{e-'^e~^ + A^e^'^ + e"^ + 7?e^) ^ CA^e^l"^, (7.9) i.^:im 

by (6.2) and (6.3). As for the particles from the family of the fugitive, note first that we have 
by Lemma 6.15, 

P(r(2 + ) = T) ^ P(T(2 + ) < r(°'^)) ^ CA^e^. (7.10) {eq:720> 

Furthermore, by Lemmas 6.14 and 6.15 and the fact that Z^^ = on the event {T = T^^^}, 

E[e-^zg\ T ^ Via^] ^ CV^P{T'^^+^ = T) ^ CAe^^^ (7.11) u,:72iy 

and likewise 

E[(e-^Z^^V, T ^ ^/^a^] ^ CVEF{t'^^+'^ = T) ^ CAe^^\ (7.12) ^e,:722> 
Likewise, if we define 

G = - Z«| ^ e^^ ^ e^/VlO, ^« < e^/(4e), Yt ^ e'}, 

then we have by (7.10), (7.11), Lemmas 6.14 and 6.15 and Markov's inequality, 

p(G'=, T ^ £3/4^3) ^ p^Qc^ J, ^ ^3/4^3 1 J. ^ 2^(i))p(r = r(i)) + p(r = r(2+)) 

^ CeAie-^e-"^ + Ae^l^ + e"^ + r/e^/e) + A^e^ (7.13) <e,:727> 

^ CA^e'^l\ 
by (6.2) and (6.3). Finally, defining 

G = {Zt^ e-^l\ Yt ^ e-^'V^}. 
we get by Lemmas 6.13 and 6.14 and Markov's inequality, 

P(G") ^ CAe^l'^. (7.14) <eq:737> 

Altogether, since Zt = Z^"-^'^) + + + ^T, we have G n G n G c Gbuik for large A, 
and thus, by (6.1), (7.9), (7.13) and (7.14), 

P(G^^uik> T ^ e^l^a^) ^ CA^e^'^. (7.15) i.^:im 
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The particles touching the right barrier after the breakout. Now, from the time T+( 
on, we are moving the barrier according to the function /a, which is equivalent to having the 
variable drift —fit = —fJ- — fA{t/cL'^)/o?- Note that on Gfug, the variable A is ^r+^-measurable 
and that Ti = T + a^/^. Since Zt ^ 2e^/e on Gbuik ^ G^fugj we now have for large a, by the 
hypotheses on the functions {fx), 

on Gbulk ^ Gfug : 1 1 /a 1 1 ^ \/a and A — /A((a^''^ — C)/^*^) = (7-16) {eq:Deita_f} 

We now show that on the good events defined above, with high probability there is no 
particle hitting the right barrier between the times T + and T + o^/^ and the descendants 
of the particles that hit the right barrier between T and T + are negligible. For this, we 
start afresh the notation of the tiers from the stopping line on, indicating this change of 
notation by a prime ('), i.e. for all particles u, such that jVt < {u,t) for some t, we set (7q{u) 
to be the second coordinate of ^Ar a (u, t) and define a'^ and by 

t'^{u) = inf{s ^ an : Xu{s) = a}, 
^n+i{u) = inf{s ^ Tn : Xu{s) = a-y+ {co- fi){s - Tn) - /a((s - (T + C))/a^)}. 

The stopping lines ^'J"^'^ etc. are then defined as in Section 6.1, adding /a to the definition. 
Note that we assumed there that / = 0, but we will not use the results from that section, 
such that there is no conflict. We then define 

/(O) „/(0) 



Gbarrier = f] ^^"''^ ^ " K+a^ = 

n {Z'ff ^ e\^} n {y}^^ ^ ST^e^} n {iig^ = 0}. 
Now, first note that on Gbuik ^ Gfug, we have by (6.3), 

Zt ^ 2e"^e'^, Yt ^ 2e"^/^ and \f{u, s) e ^ : s ^ T + C- (7.17) {eq:755> 
It follows from Lemmas 5.5 and 5.6, (6.14) and (7.17) that for large a, 

PfCbuik n Gfug n U (r(«'^))-) pbB[R'^''}^,] ^ Ge-^/^e"^, (7.18) u,:rm 



and 

P(Gbuik n Gfug, - Rf}^, ^ 1) = 0(1). (7.19) i.,:rm 

As for the tier 1 particles, if we set 



G' = Gbulk n Gfug n f| r("'^), 



(n,s)G,«'(0) , 



we have by Lemmas 5.5 and 5.6, 



E 



T + a 



2 ^^"-'^ ^ CAE lo'I^rla^ ^ CA{e-'^l^ + o{l)) ^ Ce^+^'\^, 



(n,s)s.«'(°' , 



(7.20) •[eq:758> 
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by (6.2). Equations (7.19) and (7.20) together with Proposition 5.2, Lemmas 5.6 and 5.5 and 
Markov's inequahty now give for large a, 

p(g' n {{Z'ff > e'^e^} u {vf^^ > a'^l'^e^} u {U^^ ^ 1})) ^ Ce^/^ (7.21) ,.^am 

Equations (7.18), (7.19) and (7.21), together with (6.2), (7.15) and (7.4) now prove that 

P(G'bulk ^ Gfug n G^ja^j.j.jgj.) ^ Ce^/^. (7.22) {eq:750} 

The particles that stay in the interval (0, a) after the breakout. On Gbuik ^ G'fug, we 
have Zt ^ ne^ for large A, since Z^^ ^^'^ ^ ee^ on Gfug and Zt-^^'^'^^ > Ke^-e3/2(e^ + l/4) 
on Gbuik by the hypothesis on Zq. It follows that 

1 ( Zt\ 

on Gbulk r\ Gfug : A = — log ( — ^ ] ^ 0. (7.23) {eq:767> 

By (7.16) and Proposition 5.2, we now have for large a, 

on Gbuik n Gfug : |E[zJ°^ | =^^^1 - '^e^l ^ £^^^6^/3, (7.24) <e,:77o> 

as well as 

on Gbuik n Gfug : | Var[zJ°^ | ^ .^tJI ^ C'^"^^^ + (7-25) {e,:77i> 

Equations (7.24) and (7.25) and the conditional Chebychev inequality now give for large a: 

P(Gbuik n Gfug, \zf^ - Ke^l > e^l'^e^jl) ^ C£-^l^e~^^ ^ Ce^'^, (7.26) u^-.im 
by (6.2). Hence, for large a, we have by (7.26) and (5.14), 

P(Gbuik n Gfug n Gf^J Ce^^/^^ (7.27) <e,:775> 

where 

The probability of G„. Equations (7.4), (7.15), (7.22) and (7.27) now give for large A and 
a, 

P ((Gbulk n Gfug n Gbarrier H Gi^t)^) ^ CA^e'>'^ + P(r > e3/4e^) ^ ^9/8/2, (7.28) <e,:776> 

by (6.2) and Lemma 6.2. Now note that on Gfug n Gbarrien the first and second points in the 
definition of Gi from the statement of Proposition 7.4 are verified for large a, and the third 
point is verified on Gbarrier ^ G^it for large a. In total, we have for large A and a, 

P(Gi) ^ l-e^/^ 

The statement for the probability of Gn now follows readily by induction, since conditioned 
on the event G„, the process started at the stopping time T„ satisfies the hypotheses of the 
proposition. 
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The Fourier transform of the barrier process. Define 



idrift 



{Ke^r'Zo - 1 + {Ke^)-\Z'^> + Z^^' + Z^^), and 



^jump 



and note that Adrift and Ajump are independent random variables. By tlie definition of Gbulk 
and tlie liypotliesis on Zq, we liave for large A, 



on Gbulk : A = log ( 1 + Adrift + Ajump 



0(e3/2) 



Note that by (7.7) and (7.11) and the hypothesis on Zq, we have 

«^~^co7o(^ + log e + c + oa(1)), 
and by (7.8) and (7.12) and the inequality (x + y + z)^ ^ 8(2;^ + + z^) 



E |^Adriftl(T==:v^a3) 



E 



[(Adrift)' 



L(T==:v^a3) 



0{e'A'). 



Note that (7.31) implies 

P (lAdriftl > e^l\ T ^ ^ea^) = 0{e^'^A^) = 0(e^/6), 



(7.29) -Ceq:778} 



(7.30) -[eq:Delta_drift} 



(7.31) {eq:Delta_drift_.2} 



(7.32) -[eq:Delta_drif t_bomi 



by (6.1). Now, since log(l + a + b) = log(l + a) + log((l + 6)/(l + a)), we have by (7.16), 
(7.28), (7.29) and (7.32), for large A and a. 



E[e 



(iA/co) log(l+ Adrift) - 



(iA/co)log(l+ 



(7.33) •(eq:782} 



for any A e R. We will first study the term concerning Ajump- Write Z = Z^"^^''^'^ and let p be 
a real- valued constant with \p\ < e^l^. Then, 



E[e 



where 



(iA/co)log(l+4H2P) 



] = E[e(*^/^'')'°s(^+''^^) I Z > ee^\ + O(P'^(t(0^) > C)) 
J 5(x)P(e^^Ze dx|Z > ee^) + 0(e2), 



(7.34) {eq: 



785} 



g{x) = exp I — log 1 + — — 
V Co V 1 + P 



By definition of pB and t], we have for large A and a, 

J" xP(e-^Z e dx I Z > ee^) = {pb + 0{v))-'B[Zlf^,,A^z^,,A)] 

= co7o(-loge + logK + oa(1)), 
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by (4.6), (6.2), (6.3), (6.14) and (6.13). It follows that 



r 



g{x)P{e-^Z edx\Z > se^) = 1 + a^-p^(- log e + log k + 0^(1)) 



r 



1^ So 



g{x) - 1 - iAY-p^xl(^.^«,)P(e"^Z edx\Z > ee^). (7.35) m:786} 



Now define h{x) = g{x) — 1 — zX^k^^^q/^I + p))xl(^j^^^^ for x ^ and denote by h^{x) its 
left-hand derivative. Note that \h{x)\ ^ C(l a x^) and \h^{x)\ ^ C{x^^ a x^) for x ^ 0. Now, 
by integration by parts, (4.5) and (6.13), we have for large a. 



I 



■00 

h{x)P{e-^Z e dx \ Z > ee^) 



(7.36) -Ceq :786a} 



= h{e) + p^^l + oa(1)) (J /i" (a;)P(^ > xe^) dx + {h{l) - /i(l-))P(Z > e' 

= C070 ( J /i-(x)^ dx + (/i(l) - h{l-)) + oa{1)) 

= co7o(l + Oyi(l)) h{x)^dx. 

JO ^ 

Now, one readily sees that 

r°° 1 /iA r°° \ 

/i(x)^dx = (c +oa(1) + 0(p)) + e*^"^ - 1 - iAxl(^^i) A(dx) , (7.37) {eq:786b> 

Jo ^ Jo "~ ^ 

where A(dx) is as in the statement of Proposition 7.4 and c' is a constant depending only on 
cq. Equations (7.34), (7.35), (7.36) and (7.37) and the Taylor expansion of e^^' at x = now 
yield 

j,^g(iA/co)log(l+%HP)-| 

= exp (^iA(- log e + log K + c' + OA (1) + Co J e*^^ - 1 - iAxl(^.^i) A(dx)^ +0{e\loge\p) . 

(7.38) •teq:789> 

Coming back to (7.33), we have by the Taylor expansion of (1 + x)*'^/'^^ at x = 0, 
E[e(.A/co)(log(l+A,,,)+0(.|log.|A,,a))i^^^^^3^ |A,,al<eV3)lG,.,J 

= E[(l + ^Adrift + 0(ALft))l(T^via«)] + Oie^/^) by (7.28) and (7.32) 

= 1 + iXn-^-foiA + loge + c + oa(1)) + P(r > ^/ea^) + 0{e^/^) by (7.30) and (7.31) 
= expiAK"So(^ + loge + c + oa(1) + 0(e^/®)), 

where the last equation follows from Lemma 6.2 and the Taylor expansion of at x = 0. 
This equation, together with (7.33) and (7.38) and the fact that Ajump is independent from 
Adrift) T and Gbulkj yields (7.3) in the case n = 1. For general n, we note that 
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Now, by (7.3) in the case n = 1, we have 



E[e 



iXXr 



-1 1g Je^''''''"'('^"(-^)+*'^^+°'*(^)) + 0(e^+'') 



and 



|E[e*^^-"lG„_AGJI^P(Gn-l\Gn)- 

The statement now follows easily by induction over n and the previously established fact that 
P(Gn) ^ l-nO(ei+^). 

7.2 Proof of Theorems 7.2 and 7.3 

We set 7 = 7o/(kco) and define the process {X")t^Q by 



Proposition 7.6. T/ie process {X")t^Q converges in law (with respect to Skorokhod's topology) 
to the Levy process {Lt)t^o defined in Theorem 7.2. 

Proof. Denote by (^/')t>0 the natural filtration of the process X^' , and note that = 



^"u _i I ^T. 1 • In order to show convergence of the finite-dimensional distributions, it 
is enough to show (see Proposition 3.1 in [47] or Lemma 8.1 in [32], p. 225), that for every 
A e R and t,s ^ 0, 



E 



0, 



(7.39) 



-[prop :Xt_second} 



{eq:910} 



as A and a go to infinity. Now, define n := [tj ^\ and m := [{t + s)j Then we have by 
Proposition 7.4, 

= exp ((m - n)-f{K^{A) + iXA + oa{1) + 0{e^)) - iAAs)lG„. 
= exp [s{K^{A) + oa{1) + 0(e^)))lG„, 

(7.40) {eq:912} 

because we have |(m — n)j — s| ^ 7 = A~^oa{1) + o(l), by (6.1) and (6.14). In total, we get 
for A and a large enough. 



E 



By Proposition 7.4, this goes to as yl and a go to infinity, which proves (7.39). 

In order to show tightness in Skorokhod's topology, we use Aldous' famous criterion [4] 
(see also [12], Theorem 16.10): If for every M > 0, every family of (^/')-stopping times 
r = t{A, a) taking only finitely many values, all of which in [0, M] and every h = h{A, a) ^ 
with h{A, a) ^ as A and a go to infinity, we have 



X^^^j — X" 0, in probability as A and a go to infinity, 



(7.41) 



{eq:920} 
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then tightness follows for the processes X'l (note that the second point in the criterion, namely 
tightness of X'l for every fixed follows from the convergence in finite-dimensional distribu- 
tions proved above). Now let r be such a stopping time and let Vr be the (finite) set of values 
it takes. We first note that since Gn ^ Gn+\ for every n e N, we have for every t e Vj- and 
every A and a large enough, 

P(q,^-,j)^P(G[,,^-ij) = 0(Me^). (7.42) <e,:«..> 

by Proposition 7.4. Moreover, since =^T^j^_ij for every t ^ 0, we have for every A > 0, 



teVr 



2 E[E[e^^(^"+'>-^") I + 0{Me') by (7.42) 

1 - 0{Me^)) + 0{Me^), by (7.40), 



which converges to 1 as ^ and a go to infinity. This implies (7.41) and therefore proves 
tightness in Skorokhod's topology, since M was arbitrary. Together with the convergence in 
finite-dimensional distributions proved above, the lemma follows. □ 



A coupling with a Poisson process. Let (V„),i^o be a sequence of independent expo- 
nentially distributed random variable with parameter pbc'^ttk,. In order to prove convergence 
of the processes and Xt, we are going to couple the BBM with the sequence (Vn) in the 
following way: Suppose we have constructed the BBM until time T„_i. Now, on the event 
Gn-i, by Lemma 6.6, the strong Markov property of BBM and the transfer theorem ([43], 
Theorem 5.10), we can construct the BBM up to time T„ such that 

P(Gcoupl,n) ^ 1 - 0(e2), (7.43) ie,:.2» 



where 



Geoupi,n = {|(Tr - Tn-i)/a^ " K| ^ 8^^^} n {|(rr _ r„_i)/(aVn) - 1| ^ V^) ^ e^}, 

where is here the time of the first breakout after T„_i. On the event G'^_i, we simply 
let the BBM evolve independently of (V^)j>n- Now, define 

~ ^ Gcouplj - 

Then, on G'^, we have T„ = + a^/^, whence for large A and a, 

on G'n : |(T„ - Tn-i)/a^ - K| ^ 2e^/^ and |(r„ - Tn-i)/{a^Vn) - 1| ^ (7.44) u,:.soy 

Furthermore, by (7.43) and Proposition 7.4, there is a 5 > 0, such that for large A and a, 

P{G'J > 1 - nO(e^ + '^) (7.45) {eq:prob.Gnprime> 
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Proof of Theorem 7.3. Let d denote the Skorokhod metric on -D([0, oo)) (see [32], Section 3.5). 
Let ^ be the space of strictly increasing, continuous, maps of [0, oo) onto itself. Let x,xi,X2, ■ ■ ■ 
be elements of -D([0,oo)). Then ([32], Proposition 3.5.3), d{xn,x) ^ as n — > cx) if and only 
if for every M > there exists a sequence {(fn) in such that 



and 



sup \ipn{t) -t\ ^0, 
te[0,M] 



sup \Xn{^Pn{t)) - X{t)\ 0. 
te[0,M] 



(7.46) -Ceq:940} 



(7.47) -Ceq:941} 



If (rc^)neN is another sequence of functions in D([0,aD)), with d{x'^,x) 0, then by the 
triangle inequality and the fact that $ is stable under the operations of inverse and convolution, 
we have d{xn,x) — »• if and only if there exists a sequence ((/^n) in such that (7.46) holds 
and 

sup |x„(v^„(t)) -x;(t)| ^0. (7.48) 

ie[0,M] 

For every A and a, we define the (random) map fA,a e *^ by 



-Ceq:942> 



TT 



ipA,ait) = ((1 - r)Tn + rTn+i)^, if t = 7(n + r), with n e N, r e [0, 1]. 



Let M > and define tt-m = [-^t]- Then we have 



sup \^pA,a{t) — t\ ^ max 
te[0,M] ' ne{0,...,nM} 



TT 



Tn - in 



and 



sup \X'-I — X' \ ^ max A 

te[0,M] ^^'"^^ ne{0,...,n,,} 



TT 



Tn - in 



Now note that 7 = ^£^[Vi], and by Doob's L?' inequality we get 



(7.49) {eq:944} 



(7.50) -Ceq:945> 



PI max 

V ne{0,...,njv/} 



i=l 



> ^ 4e-2/3nA/ Var(V-) = 0{e^'^). 



Furthermore, on the set G^^^^, we have for every n ^ tt-M; 



«=0 



In total, we get with (7.49) and (7.50), as A and a go to infinity. 



VM > : sup \ipA,a{t) - t| V \X'I - X' 



te[Q,M] 



'PA,a{t) ' 



0, in probability, 



which is equivalent to 



2 2-^4l A ( sup \ipA,a{t)-t\ V 



(7.51) 



-Ceq:948} 



0, in probability. (7.52) {eq:949> 
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Now, suppose that A and a go to infinity along a sequence (A„, an)neN and denote by X'^^ 
^'L a„ and 

fA„,a„ the processes corresponding to these parameters. By Proposition 7.6, and 
Skorokhod's representation theorem ([12], Theorem 6.7), there exists a probabihty space, on 
which the sequence [X'j^^ converges almost surely as n — > OD to the limiting Levy process 
L = [Lt)t^Q stated in the theorem. Applying again the representation theorem as well as 
the transfer theorem, we can transfer the processes X'j^^ and ^PA„,a„ to this probability 
space in such a way that the convergence in (7.52) holds almost surely, which implies that 
the convergence in (7.51) holds almost surely as well. By the remarks at the beginning of the 
proof, it follows that on this new probability space, 

«„,a„,^) ^ d(X;,„,.„,X^„,.J + d{X'l^,^,L) - 0, 
almost surely, as n ^ cc. This proves the theorem. □ 

Proof of Theorem 7.2. By the virtue of Theorem 7.3, it suffices to show that for every < 
ti < t2 < ■ ■ ■ < tk we have 

P(Vi : X[^J = 4,3) ^ 1. (7.53) <e,:950> 

Let n := \2{tk + 2)/E[ViW, such that n = 0{e^^), by (6.14). By Chebychev's inequality, we 
then have 

n n 

P(2 V- ^ tfc + 2) ^ P(2(^* - ^[^']) ^ = 0(nVar(y,)) = Oie). (7.54) ^e,:952> 

i=l i=l 

Furthermore, define the intervals li = ti + [— 2ne^''^ — a^^^"^, 2ne^^'^], i = 1, . . . ,k and denote by 
^ the point process on the real line with points at the positions Vi , Vi + V2 , Vi + V2 + V3 , . . . . 
Then ^ is a Poisson process with intensity l/£'[Vi] = 0{e^^) and thus, 

k 

P(=!^ n y ^ 0) = 0(e^/2^ + 0(1). (7.55) {eq:954> 

2=1 

We now have 



P(yi : Xl^J = j ^ P(t{i,j) : ha^ e [T, - rf_0, T,] j by definition 

n k 

^ Pfc;, Y^^^>tk + 2, ^ n (J = 0) by definition of G'^ 
^ 1 - O(e^) - 0(1) by (7.45), (7.54), (7.55). 

Letting A and a go to infinity yields (7.53) and thus proves the theorem. □ 
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